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The  third-order  elastic  constants  C ^ ^ j of  four  types  of  fused  silica 
were  measured  as  a function  of  temperature  between  room  temperature  and 
3°  K.  The  room  temperature  values  were  determined  absolutely  by  measuring 
the  harmonic  distortion  of  an  initially  sinusoidal  ultrasonic  wave 
propagating  through  the  sample.  Relative  values  were  measured  at  lower 
temperatures.  The  results  indicate  that  a relatively  large  OH  content  in 
fused  silica  gives  rise  to  a relatively  large  value  of  an<*  **,at 

is  a weak  function  of  temperature  regardless  of  OH  content. 

The  relationship  between  the  experimentally  defined  ultrasonic 

*^11  + ^111 

nonlinearity  parameter  6 = - and  the  Grtlneisen  parameter  y 

also  has  been  determined. 
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CHAPTER  I 


INTRODUCTION 

A number  of  the  physical  properties  of  solids  result  from  the  fact 
that  solids  are  inherently  nonlinear.  Among  the  properties  are  thermal 
expansion,  attenuation  of  high  frequency  sound  waves,  heat  conduction, 
pressure  variations  of  second-order  clastic  constants,  and  waveform 
distortion  of  sound  waves  passing  through  a solid.  These  nonlinear 
effects  arise  because  of  anharmonicity  of  the  interatomic  potential 
(the  deviation  from  a potential  that  is  quadratic  in  the  interatomic 
displacement  from  equilibrium). 

In  the  present  experiments  anharmonicity  is  investigated  from  the 
distortion  of  ultrasonic  waves  passing  through  a solid,  and  the  results 
are  related  to  the  Griineisen  parameter  which  is  a thermodynamic  measure 
of  nonlinearity  in  solids. 

Previously  the  investigation  of  anharmonicity  in  solids  has  been 
pursued  from  the  more  purely  thermodynamic  investigations.  Thermal 
expansion  has  been  among  the  most  measured  of  the  thermodynamic 
properties,  and  data  for  many  materials  have  been  available  for  many 
years.  For  example,  Fizeau  (1863)  was  among  the  earliest  investigators 
to  measure  thermal  expansion. 

Griineisen  (1926)  developed  the  relationship  between  the  thermal 
expansion  coefficient  and  a parameter  which  was  assumed  to  be  temperature 
independent  (the  Griineisen  parameter  y)  . Barron  (195S)  and  Sheard  (1938) 
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wore  among  the  first  to  make  calculations  relating  the  GrVinciscn 
parameter  to  clastic  data.  They  limited  their  calculations  to  the  cases 
of  high  and  low  temperatures.  Collins  (1963),  Schuele  and  Smith  (1964), 
Brugger  and  Fritz  (1967),  and  C.erlich  (1975)  made  more  extensive  calcu- 
lations to  include  the  temperature  dependence  of  the  Gruneisen  parameter 
from  elastic  data,  using  the  quasiharmonic  Debye  model  of  solids.  Thus, 
the  Gruneisen  constant  is  a convenient  link  between  thermal  expansion 
data  and  elastic  data. 

The  agreement  of  the  calculated  dependence  of  the  Gruneisen 
parameter  as  a function  of  temperature  with  values  measured  from  thermal 
expansion  has  been  limited,  because  of  the  lack  of  availability  of 
elastic  data  (particularly  third-order  elastic  constants)  as  a function 
of  temperature.  For  many  materials  such  as  Cu,  Ag,  Au,  Al,  and  Na 
qualitative  agreement  exists  hut  some  materials  exhibit  an  anomalous 
behavior  of  the  thermal  expansion  coefficient  as  a function  of  tempera- 
ture. One  such  material  is  fused  silica.  Fused  silica  has  a positive 
thermal  expansion  coefficient  until  it  reaches  a temperature  between  .3 
and  .4  of  its  Debye  temperature.  In  that  range  the  thermal  expansion 
coefficient  becomes  negative  and  continues  to  grow  more  negative  as  the 
temperature  is  lowered  (White,  1964).  Several  investigations  (Clark 
and  Strakna,  1962;  White,  1964;  Klein  and  Mountain,  1962;  McCammon  and 
White,  1963)  indicate  that  at  least  part  of  the  explanation  for  this 
behavior  of  fused  silica  may  be  due  to  the  dominance  of  transverse 
acoustic  modes  at  these  temperatures;  however,  this  hypothesis  can  be 
tested  only  as  information  on  the  behavior  of  third-order  elastic  (TOE) 
constants  as  functions  of  temperature  is  collected.  Since  the 
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third-order  clastic  constants  determine  how  a large  amplitude  ultrasonic 
wave  distorts  as  it  propagates  through  the  solid,  TOE  constants  can  be 
calculated  from  ultrasonic  wave  distortion. 

The  distortion  of  an  ultrasonic  wave  as  it  passes  through  a solid 
may  be  studied  by  use  of  either  classical  continuum  mechanics  or 
quantum  mechanics.  Classically,  this  phenomenon  may  be  described  by  a 
nonlinear  wave  equation  whose  solution  predicts  the  growth  of  higher 
harmonics  as  the  wave  propagates  through  the  solid.  Quantum 
mechanically,  this  is  equivalent  to  two  phonons  of  one  frequency  inter- 
acting to  produce  one  phonon  at  the  second  harmonic  of  that  frequency. 
Both  descriptions  depend  upon  the  introduction  of  terms  of  order  higher 
than  quadratic  in  the  strain  potential  energy  function  of  the  solid. 

The  coefficients  (elastic  constants)  of  the  power  series  expansion  of 
the  strain  potential  energy  function  are  related  to  the  coefficient 
(force  constants)  of  the  power  series  expansion  of  the  interatomic 
potential  function.  For  example,  Coldwel 1-Horsfal 1 (1963)  gave  a 
relation  between  force  constants  between  molecules  in  a cubic  crystal 
and  the  elastic  constants  of  the  solid,  and  Ghate  (1965)  used  the  Born 
model  of  ionic  solids  to  compute  the  third-order  elastic  constants  of 
alkali  halides. 

Several  experimental  techniques  have  been  used  to  determine  the 
nonlinear  behavior  of  solids.  Hughes  and  Kelley  (1953)  were  first  to 
obtain  the  complete  set  of  third-order  elastic  (TOE)  constants  for 
isotropic  materials  by  measuring  the  variation  of  ultrasonic  wave 
velocity  with  hydrostatic  pressure  and  with  uniaxial  stress  in  poly- 
styrene, iron,  and  pyrex  glass.  Skove  and  Powell  (1967)  pointed  out 


5 


measurement  of  the  absolute  amplitude  of  the  fundamental  and  harmonics  of 
the  ultrasonic  waveform  after  propagation  through  the  solid  and  calculated 
the  TOli  constants  of  Copper.  Yost  and  Breazeale  (1973)  used  this 
receiver  to  measure  fused  silica  and  combined  the  results  with  those  of 
Dunham  and  Huntington  (1970)  to  obtain  the  first  complete  set  of  truly 
adiabatic  TOO  constants. 

Mackey  and  Arnold  (1969)  and  Meeks  and  Arnold  (1970)  used  this 
method  to  measure  combinations  of  TOf:  constants  from  room  temperature 
to  100  °K.  Peters,  Breazeale,  and  Par6  (1968)  developed  a pneumatically 
controlled  variable  gap  capacitive  receiver  capable  of  making  accurate 
measurements  at  lower  temperatures.  Peters,  Breazeale,  and  Par£  (1970) 
measured  combinations  of  TOE',  constants  of  copper  to  liquid  nitrogen 
temperatures,  and  Yost  and  Breazeale  (1974)  made  similar  measurements 
on  germanium.  Bains  and  Breazeale  (1976)  extended  the  capability  of 
this  technique  to  liquid  helium  temperatures  for  measurements  on 
germanium. 

In  the  experiments  described  in  this  dissertation  the  apparatus 
of  Bains  (1974)  was  used  to  measure  harmonic  distortion  of  a 30  MHz 
ultrasonic  wave  in  four  types  of  fused  silica.  The  TOP  constants  arc 
measured  between  room  temperature  and  3 °K,  and  their  behavior  as  u 
function  of  temperature  is  related  to  that  of  the  Brunei  son  parameter. 

In  Chapter  II  the  theory  of  ultrasonic  nonlinearity  in  solids  is 


discussed,  and  the  relationship  between  measured  quantities  and  T01 
constants  is  derived.  The  relationship  between  the  TOP  constants  and 
the  Gruneisen  parameter  is  examined  in  Chapter  III.  Apparatus, 


procedure,  and  samples  are  discussed  in  Chapter  1\ . I’inall) , Chapter  \ 
presents  results  of  measurements  on  fused  silica,  and  discussion  of  the 
implications  of  these  measurements. 
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C1IAPTIU  II 

Til  LORY  OP  ULTRASONIC:  NONI.INP.ARITY  IN  SOLIDS 

The  theory  of  a finite  amplitude  ultrasonic  wave  propagating 
through  a solid  medium  may  be  approached  either  from  a quantum  mechanical 
point  of  view  or  from  classical  continuum  mechanics.  Quantum  mechanically 
we  may  consider  three-phonon  interaction  processes  to  calculate  transition 
probabilities  using  appropriate  interaction  Hamiltonians.  These  transition 
probabilities  yield  the  corresponding  classical  results  when  integration 
is  performed  over  the  appropriate  interaction  volume  (Bajak,  1P7M. 

The  interaction  volume  is  typically  more  than  twenty  orders  of 
magnitude  greater  than  primitive  cell  (intcrmolecular)  volumes  and  for 
this  reason  the  approach  of  continuum  mechanics  is  justified.  Hence, 
the  distortion  of  an  ultrasonic  wave  as  it  propagates  through  the  solid 
can  he  described  by  a solution  to  the  nonlinear  wave  equation  obtained 
from  Lagrange's  equation  for  a continuum.  The  nonlinear  terms  in  the 
wave  equation  enter  in  two  ways.  First,  the  power  series  expansion  of 
the  potential  energy  density  (strain  energy  density)  in  the  Lagrangian 
contains  terms  through  the  third  order.  These  third  order  terms  lead 
to  nonlinear  terms  in  the  wave  equation.  Second,  the  strains  contain 
nonlinear  terms  themselves  since  they  represent  a finite  deformation  of 
the  propagation  medium.  Solution  of  the  nonlinear  wave  equation  yields 
relationships  between  the  harmonic  amplitudes  of  the  distorted  ultrasonic 
wave  and  the  coefficients  of  the  strain  energy  density  expansion  (the 


elastic  constants). 
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I.  1)1:1*  IN  I TION  OF  LACRANCIAN  STRAINS 


A point  in  an  unstrained  solid  may  he  represented  by  the  vector  a. 
lVher»  the  solid  is  strained  by  application  of  some  stress  the  point  is 
located  by  a new  vector  x.  The  point  thus  undergoes  a displacement 

u = x - a (II.  1) 


when  the  stress  is  applied.  The  vector  a describes  the  point  in 
Lagrangian  (laboratory)  coordinates  and  the  vector  x describes  the  same 
point  in  Eulerian  coordinates  (which  arc  fixed  in  the  solid).  The  strain 
associated  with  the  deformation  of  the  solid  may  be  described  either  in 
Lagrangian  or  Eulerian  coordinates.  Since  we  must  make  experimental 
measurements  in  a laboratory  reference  frame  we  shall  use  Lagrangian 
coordinates. 

Consider  two  points  in  the  solid  in  the  unstrained  state  and  lot 
the  distance  d£  between  them  be  given  by 


dfc“  = | da | “ = daT  (11.2) 

where  Einstein  notation  (repeated  indices  are  summed)  is  assumed.  The 
distance  between  the  same  two  points  in  the  strained  state  dfi.*  is  given 
by 

dl'2  = |dx  | 2 = dxT  . (II. 3) 

From  Eq.  (It.l)  we  may  write 

2 2 ^Ui  2 

dP, * = (da.  + du.)  = (da.  + - — da.) 

ii  l aa . i 

J 


(II. 4) 


■p 


-4 


Irom  Pqs.  (11.2)  and  (ir.4), 


2 *> 

dsi * ^ - dK4*  = 2n-  .da. da. 

ij  » J 


( ri .5) 


where  n—  arc  the  elements  of  the  Lagrangian  strain  tensor  defined  by 


ij  2 


I 

9u. 

l 

2 

9a . 
) 

1 Z* 

i.  9a., 
i 


(II. 6) 


lor  small  amplitude  strains  the  second  order  terms  in  Eq.  (11.6)  are 

j/-3u.  3u.' 

commonly  dropped  and  the  n.  • reduce  to  the  form  -=■  - — - ♦ — used  in 

ii  2 ( aa  . 3a . 

) i 

linear  theory  (Green,  1973).  However,  they  are  essential  in  the  proper 
description  of  the  subsequent  nonlinear  effects. 

IT.  TUP.  EQUATION  OF  MOTION 

The  equation  of  motion  in  laboratory  coordinates  can  be  obtained  from 
Lagrange's  equation  for  a continuous  medium  (Goldstein,  1968), 


dl 

9L  1 

4. 

’ d 

9L 

Lik_ 

dti 

l3*jj 

K 

9x . 

**  r 1a 

j 9x . 
; 1 

= 0 


(II. 7) 


9a, 


where  the  Lagrangian  density  l is  given  by 


1 . 2 

L = 2 P0*.  - Mn)  , 


(II. 8) 


Pq  is  the  mass  density  of  the  unstrained  medium,  and  4>(n)  is  the  strain 
energy  density.  It  is  seen  that  the  strain  (and  hence  the  nonlinearity) 
enters  the  equation  of  motion  only  through  the  strain  energy  density. 
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II J.  Till:  STRAIN  HNHRC.V  DHNSITY  AND  FUSTIC  CONSTANTS 


When  an  clastic  body  is  deformed  the  agent  applying  the  stress 
necessarily  docs  work  on  that  body.  If  we  combine  the  first  and  second 
laws  of  thermodynamics  we  can  write  the  change  in  internal  energy  of 
the  deformed  body  as 


dU  = TdS  - dW 


(II. 9) 


where  U is  the  internal  energy  function  per  unit  volume,  T is  the 
temperature,  S is  the  entropy  and  dW  is  the  work  done  by  the  agent  in 
deforming  the  body.  Assuming  that  the  applied  stress  on  the  body  is 
represented  by  the  stress  tensor  Thurston  (1964)  showed  that 
licj.  (1 1. 9)  can  be  written  as 


dU  = TdS  + — — t . . dn. • . 

P0  U ’J 


(II. 10) 


Thus,  we  may  consider  the  internal  energy  function  IJ  to  be  a function  of 
the  entropy  S and  Lagrangian  strains  n— . 

If  we  expand  U(S,n^j)  in  a power  series  of  the  strains  we  get 


2 

U(S,n.J  = U(S)4P  -11-  n..n.t 
U 2 °piij3n].Js  ij  ki 


31  0 3nij3rWv)s  ij  k* mn 


(ii. m 


Writing  the  Helmholtz  free  energy  per  unit  volume  F = U - TS,  the 
differential  of  F is 


«***  > * ' - ■!?*»&*% 


dF  = -SdT  •*■—!.  ,dn.  . . 
P0  U 
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(11.121 


llcr.ce,  we  may  consider  F to  be  a function  of  T and  n^j,  and  expand  F in 
a power  series  analogous  to  Eq.  (11.11), 


•? 

3"F 


l-TT.n  ) = >-(T)  * f P()| 

l ij  Ml 


nijnU 


JT  po 


3 

3 F 


. 3n • . 9n»  „ 9n 
( ij  *k£  mn 


ni  i nk£r'mn 


(11.13) 


Brugger  (1%4),  making  use  of  the  strain  dependence  of  the 
thermodynamic  state  functions  II  and  F,  defined  the  coefficients  of  the 
nth  order  in  the  power  series  in  Eqs.  (1 1. 11)  and  (11.13)  to  be  the 
adiabatic  and  isothermal  clastic  constants  of  nth  order,  respectively. 
Thus , 

rS 

C i jk£  ...  P0 

and 

rT 

ljki  . . . p0 

If  we  equate  the  strain  energy  density  in  Oq.  (1 1.8)  to  cither  the 
internal  energy  per  unit  volume  IJ  or  the  Helmholtz  free  energy  per  unit 
volume  F , respectively,  then  the  equation  of  motion  (Fq.  (11.7) | that 
results  would  describe  adiabatic  or  isothermal  notion,  respectively. 


3nij3nu 


BniJ  fadiabatic 

' ( n “ 2 . 


(11.14) 
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In  general,  there  are  81  second-order  and  729  third-order  elastic 
constants.  Lattice  symmetries  of  the  medium  can  reduce  the  number  of 
independent  constants  so  that  in  isotropic  media  there  are  only  two 
independent  second-order  elastic  constants  and  three  third-order  elastic 
constants  (Birch,  1947), 

Since  the  Lagrangian  strains  defined  by  F.q.  (II. 6)  are  symmetric, 
only  six  of  each  set  of  nine  strains  are  independent.  Hence,  we  may 
introduce  the  following  contraction  of  subscript  notation  (Voigt,  1928): 

11  ->  1 , 22  » 2,  33  -*  3,  23  - 4,  13  5,  12  -**  0 . 

If  we  combine  Bruggor's  definitions  of  the  elastic  constants,  Voigt’s 
contracted  notation,  and  Birch's  study  of  crystalline  symmetries,  then  the 
two  independent  second-order  elastic  constants  for  isotropic  solids  may  be 
written  as 


12 


‘44 


w i t h 
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11 
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12 


2C 


44 


and  the  three  independent  third-order  constants  as 


r r r 
111*  144'  456  * 

IV.  THE  ADIABATIC  NONLINEAk  WAVE  EQUATION 

If  we  let  the  strain  energy  density  in  Eq.  (11.8)  he  the  internal 
energy  function  per  unit  volume  I given  by  Eq.  (11.11)]  and  substitute 
Eq.  (If. 8)  into  Lagrange's  Eq.  (1 1.7)  we  get 


Equation  (11.19)  is  the  adiabatic  equation  of  motion  (and  Cjj  and  Cjj^ 
arc  adiabatic  elastic  constants)  since  the  internal  energy  function  was 
used  for  the  strain  energy  density,  For  the  frequency  range  used  in 
these  experiments  the  use  of  the  internal  energy  function  is  justified. 
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V.  SOLUTION  TO  THE  NONLINEAR  WAVE  EQUATION 


Breazeaie  and  Ford  (1965)  solved  the  equation  of  motion 
|Hq.  (11.19))  by  assuming  an  initial  sinusoidal  disturbance  /.sinwt  at 
a = 0.  Their  solution  may  be  written  for  isotropic  solids  as 


1 1 + 2 2 

U = A_sin(ka  - tot)  - ( 57= -)  A.  cos  2(ka  - wt) 

1 8Cn  1 


(11.20) 


where 


A is  the  amplitude  of  the  fundamental, 
k is  the  propagation  constant,  2tt />,  and 
a is  the  sample  length. 

Inspection  of  Hq.  (11.20)  reveals  that  the  initially  pure  sinusoidal 
disturbance  becomes  distorted  as  it  propagates  through  the  medium.  That 
is,  a wave  of  frequency  2io  is  generated  and  grows  linearly  with  propa- 
gation distance  in  the  sample.  It  is  interesting  to  note  that  if  one 
considers  a ’’linear”  solid  to  be  one  in  which  no  second  harmonic  is 
generated,  then  the  necessary  condition  is  that 


111 


= -3C 


11 


(II. 21) 


and  not  that  the  third-order  elastic  constants  be  zero.  So  far,  no 
material  has  been  found  which  satisfies  the  condition  of  Eq.  (11.21). 


VI.  THE  NONLINEARITY  PARAMETER 


In  Eq.  (11.20)  me  quantities  that  arc  measured  in  t he  laboratory 
are  the  amplitudes  of  the  fundamental  and  second  harmonic  at  the  end  of 


1 5 

the  sample  opposite  that  at  which  the  initial  sinusoidal  disturbance  is 
generated.  Hence,  the  coefficients  of  the  trigonometric  functions  in 
Kq.  (11.20)  yield  a value  A^  for  the  fundamental  amplitude  and 


3C 


A,  = -(■ 
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111 
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2 2 

A^  k a 


(11.22) 


for  the  second  harmonic  amplitude.  Solving  for  we  get 


111 


2 1 

2 2 
A k a 


1) 


(11.23) 


All  quantities  on  the  right  hand  side  of  this  equation  can  he  measured, 
and  C.jj  can  be  determined. 

The  quantity 


8^ 

3 


K“a 


(11.24) 


in  Kq.  (IT. 23)  is  a significant  quantity  in  this  study.  If  6 were  zero, 
then  we  see  that  the  condition  of  Kq.  (11.21)  is  satisfied  for  no  second 
harmonic  generation.  As  we  shall  see  (Chapter  III),  8 is  related  to  the 
Gruneiscn  parameter  which  is  a thermodynamic  measure  of  material 
anharmoni city . 

We  shall  refer  to  8 as  the  ultrasonic  nonlinearity  parameter.  For 
isotropic  solids  8 is  expressed  in  terms  of  the  elastic  constants  as 


3Cn  + cm 


3C 


11 


6 * -( 


(11.25) 


chapter  m 

THE  GRUNEISHN  PARAMETER  AND  ITS  RELATION  TO 
THE  NONLINEARITY  PARAMETER 


The  parameter  $ defined  in  Eq.  (11.24)  offers  a direct  measure  of 
material  anharmonici ty  through  measurement  of  elastic  nonlinearity  of 
the  material  structure.  However,  anharmonicity  can  he  investigated  from 
other  measurements.  Among  them  are  isothermal  compressibility  variations 
as  a function  of  pressure,  the  thermal  expansivity  as  a function  of 
temperature,  and  variations  of  refractive  index  with  temperature  (Primak, 
1975).  Among  those  investigating  anharmonicity  from  these  techniques 
it  has  become  common  to  express  the  results  in  terms  of  t lie  Grtlneisen 
parameter  y,  defined  by 


a « 

kt  cv  = Ks  CP 


(iir.i) 


where  a is  the  thermal  volume  expansivity,  and  the  isothermal  and 
isentropic  compressibilities,  respectively,  and  L ^ and  Cp  the  isochoric 
and  isobaric  heat  capacities,  respectively  (Grtlneisen,  1926).  It  is 
the  purpose  of  this  chapter  to  show  the  relation  of  the  Griineisen 
parameter  to  the  nonlinearity  parameter  3. 


I.  RELATION  BETWEEN  THE  STATISTICAL  MECHANICAL  MODEL  OE  A 
SOLID  AND  THE  THERMODYNAMIC  STATE  FUNCTIONS 


A solid  can  be  represented  as  the  sum  of  two  entities:  an  array 

of  atoms  at  rest  in  their  mean  positions  with  static  potential  energy 
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W 


- S..,  . 


plus  an  assembly  of  lattice  vibrations  represented  as  a set  of  weakly 
interacting  but  distinguishable  simple  harmonic  oscillators.  The  set 


of  simple  harmonic  oscillators  may  be  designated  by  a branch  index  p 
wave  vector  q,  frequency  w(p, <|)/2?r,  and  energies 


n(p,q)  = [n(^  * 


(HI. 2) 


where  n(p,q)  = 0,  1,  2,  3,  ...  and  h is  Planck's  constant  divided  by 
27 r.  The  partition  function  for  such  a system  is  (Slater,  1939) 


n(p,q) 

kT 


(HI. 3) 


pq  n(p,q) 


where  k is  the  Boltzmann  constant  and  T is  the  temperature . 

The  relation  between  the  statistical  mechanical  model  of  a solid 
and  th'  thermodynamic  state  functions  U (internal  energy)  and  F (Helmholtz 
free  energy)  is  given  by  (Hill,  1962) 


IT  k i ,32,  _ , T2 ,3£nZ, 

u = — c3T}v = kT  (1T]V 


(III. 4) 


F = -kTtnZ  . 


(1M.5) 


IT.  THE  "STRAIN"  GENERALIZED  GRIJNEISEN  PARAMETERS 


The  thermal  expansivity  a given  in  Eq.  (III.  1 ) may  be  generalized 


to  tensorial  form  by  (Thurston,  1964) 


ajk  " 3T 
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whore  n..  arc  components  of  the  Lagrangian  strain  tensor,  T is  t !te 
J k 

temperature,  and  the  subscript  i means  the  thermal  strain  derivative 
is  evaluated  at  constant  stress.  In  terms  of  the  Helmholtz  free  energy 
we  nay  write  (Brugger  and  Fritz,  1967) 


a..  = - [ S! 
3k  ^ 3 


.T  3T/V 
5ikrs  3T3n 


(111.7) 


where  the  Sjkrs  are  the  isothermal  elastic  compliance  coefficients  and  V 
is  the  volume  of  the  solid. 

If  we  assume  that  $ and  w(p,q)  depend  on  temperature  only  through 
the  lattice  dimensions  (quasi harmonic  approximation)  (Liehfried  and 
Ludwig,  1961),  then  from  Fqs.  (III. 3)  and  (III. 5)  we  may  write 
liq.  (1 1 1.7)  in  the  form 


> ' k k 'W'slp'5,c,M'> 


(111.8) 


where  the  strain  generalized  Grtlneisen  parameters,  Yrs(p,q),  are 
defined  bv 


Yrs(P,q)  = 


1 

3w(p,q) 

w(p,q) 

3nrs 

(HI.9) 


These  strain  generalized  Grtlneisen  parameters  are  indirectly  dependent  on 
temperature  through  the  strains  n^s*  The  C(p,q)  are  the  modal  heat 
capacities  given  by 


«)  = 


(p,q) 

uMlVil,-’  t.  It 
_J kl  

ht.i(p,q) 

, . kT  .,2 


(I 11.10) 


They  are  directly  temperature  dependent 


III.  THE  "VOLUME"  GENERAL  1 2Hl>  CRUM:  I SEN  PARAMETERS 


The  thermal  expansivity  a may  he  expressed  in  terms  of  the  tensorial 


ri . . by 
rk  ' 


= l a., 
j JJ 


an. in 


Substituting  Eq.  (III. 8)  into  Eq.  (II 1. 11)  gives 


“ “ l (?  I S{irs 

->  vj  rs  1 


(III. 12) 


If  we  define 


T(P.q')  = ^ l £ S..rs  Yrs  (p,q 


(111.13) 


where  is  the  isothermal  compressibility  and 


Y(p,q)  = - 


V a»*»(p,q) 

, * av  t 

w(p,q) 


(III. 14) 


then  Eq.  (1 1 1. 12)  may  be  written  as 
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a = 1 Y(P»q)C(p,q)  . (111.15) 

PM 

liquation  (111.15)  is  the  form  for  a Riven  by  Slater  (1959)  with  vip,q) 
being  defined  by  liq.  (1 1 1. 14).  To  distinguish  the  y(p,q)  from  the  s t ra  l n 
generalized  Gruneisen  parameter  >rs(p,q)  of  fiq.  (II  1.9)  we  shall  call  the 
y(p,q)  of  bq . (III. 14)  the  volume  generalized  GrUncisen  parameters, 
liquation  (111.15)  then  gives  the  relation  between  the  two  parameters. 


iv.  Tur.  grUnpishn  paramftpr  and  its  relation  to  thi: 

Gf-NfRALIIflD  GRUNEISEN  PARAMETERS 
Substituting  Fiq.  (III.  15)  into  P.q.  (II  1.1)  gives 


}.  r(p,q)C(p,q) 


l C(p,q) 

P »M 

where  we  have  defined 

Cv  = l C(p,q)  . 
P>M 


(111. KO 


(II 1.17) 


Thus,  from  fiq.  (1 1 1. 16)  the  GrUncisen  parameter  y is  equal  to  a weighted 
average  of  the  volume  generalized  Grunciscn  parameter  y(p,q).  The 
weighting  functions  are  simply  the  modal  heat  capacities  C(p,q). 

Substituting  liq.  (111.15)  into  fiq.  (II  1. 16)  gives  the  dependence  of 


the  Gruneisen  parameter  on  the  strain  generalized  Gruneisen  parameter, 


, l-l  l sjjrsVs(p-‘>)Ctp^1 
_ _L  p,q  ) « •• 

KT  l C(p.q) 

p.q 
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(111.18) 


V.  CONTINUUM  MODELS  AND  RESULTING  EXPRESSIONS 
FOR  THE  GRdNFISEN  PARAMETER 


Until  now  no  particular  model  has  been  assumed  to  describe  the 
frequency  distribution  of  the  simple  harmonic  oscillators.  The  intro- 
duction of  a frequency  distribution  has  the  effect  of  assigning  relative 
"weights”  to  the  modal  heat  capacities  to  reflect  the  relative  number  of 
independent  simple  harmonic  oscillators  in  a particular  mode.  This  in 
turn  assigns  relative  "weights,”  via  the  modal  heat  capacities,  to  the 
Grtlneiscn  parameter  y. 

The  continuum  model  is  introduced  through  the  following  assumption 
(Bruggcr  and  Frit 2,  1967): 

1.  The  excitation  of  optic  modes  is  neglected  so  that  the  branch 
index  p takes  on  only  the  values  1,  2,  3 to  reflect  the 
longitudinal  and  two  transverse  acoustic  branches. 

2.  The  acoustic  modes  obey  the  Debye  frequency  distribution 
function  (per  unit  volume) 


q(P,q)dqd£l  = — —■  q2dqdtt  . 

(»r; 


nu.  i9) 


3.  The  maximum  value  of  q along  any  direction  equals  the  Debye 
radius 


% ' 


/ 7I1/3 

671' 


(III. 20) 


where  V is  the  volume  of  the  primitive  cell  of  the  solid.  The 
o 1 

volume  cf  the  Debye  sphere  equals  that  of  the  first  Brillouin 


4.  The  acoustic  modes  are  either  nondispersive  (Debye  model)  or 
they  obey  the  sinusoidal  dispersion  relation  (Born-von  Karman 
model ) 


■n  (i 

un  v — L . 

* 1 1 


(111.21 


5.  The  strain  generalized  GrUncisen  parameters  and  the  volume 

generalized  Grtlneisen  parameters  are  independent  of  wave  number 
(frequency)  in  both  models. 

The  basic  difference  between  the  dispersive  and  nondispersive  models 
is  that  the  Born-von  Karman  model  assigns  a greater  "weight"  to  the  low- 
frequency  acoustic  modes  than  docs  the  Debye  model. 

With  these  assumptions  F.q.  (111.16)  becomes  (Brugger  and  l ritr,  lDo?) 


I I dfiY(P,N)C(p,N) 

E _ 

l I dnC(P.N) 

p 
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Ccp.ni  * J d?  -y  e , 
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0 (cx  - 1) 


(111.23) 


5 = -il— 


(111.24) 


For  the  Debye  model 


Vbye^N’f- 


= ,f>  (P»N)  j 


lor  the  Bom- Von  Karman  model 


in  both  models 


0(p,N)  = -f-  S(p,N) 


(111.25) 


Vrn(P’fi’f’>  =I  (£CMi]sin(-  d . (fU.261 


(111.27) 


is  the  characteristic  temperature  of  the  mode  specified  by  p and  N,  and 


s(p.n)  = 


(III. 28) 


is  the  acoustical  wave  speed  in  the  solid. 

In  terms  of  the  modal  characteristic  temperatures  the  Debye 
temperature  0 is  given  by 


1 _ 1 c r dU 

^3  12tt  ^ ^ * .3  ' 

0 p [0(p,N) ] 


(111.29) 


VI.  THE  RELATION  OF  GENERALIZED  GRUNEISEN  PARAMETERS 
TO  ELASTIC  CONSTANTS 


Equation  (1 11.28)  implies  that  for  the  Debye  model  the  frequencies 
arc  related  for  any  state  of  strain  to  the  wave  speed  and  some  dimension 


L of  the  solid  by  ( Brugge r,  1965): 


Eq.  (111.33)  art'  evaluated  at  zero  strain,  t he  s t ra  i n genera  1 i zed 
Gruneisen  parameters  may  be  written  as 


Yrs(P.N)  « l—  [2wCp,N)U  U 

rS  2w(p,N) 


+ C . 1J  l)  )N  N 1 
rsmn  rskmunv  u v it.  n1 


(III. 34a) 


wfp  N)  = C N N IJ  l) 

1 munv  n n u v 


(111.35) 


Equation  (III. 34a)  is  the  expression  relating  the  strain  generalized 
Gruneisen  parameters  to  the  elastic  constants. 

The  volume  generalized  Gruneisen  parameters  in  terms  of  the  elastic 
constants  may  be  obtained  by  substituting  Eq.  (II 1. 34a)  into 
Eq.  (III. 13). 


y(p,N)  = - 


2KTw(p,N) 


S..  [dw(p,N)ll  D 
jjrs  1 1 9 * r s 


+ (C  + C tl  11  )N  N 

rsmn  rsmunv  ti  v m n 

(Einstein  notation). 
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VII.  THE  RELATION  OF  THE  GRUNEISEN  PARAMETERS  TO 
THE  NONLINEARITY  PARAMETER 


In  an  isotropic  solid  all  directions  are  equivalent.  Hence,  we  can 
simplify  Eq.  (111.34a)  for  the  strain  generalized  Gruneisen  parameters 
by  assuming  Nj  = l,  N2“N3=0>  and  = U7=tl5=0,  for  thc  propagation  of  a 

pure  mode  longitudinal  ultrasonic  wave  in  an  isotropic  solid. 
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Equation  (III. 34a)  then  gives  for  the  tensor  element  (r=l,  s=l)  for  al 1 
longitudinal  modes  (p^,N) 

* 3C  + Cjji 

yu(p*’n)  = - — ic-71-  un'371 

where  the  notation  is  used  to  emphasize  the  longitudinal  acoustic 
mode,  from  fiq.  (11.25)  of  Chapter  II  and  l*q.  (IJI.37),  the  relation 

A 

between  T^Cp^N)  and  the  ultrasonic  nonlinearity  parameter  3 is 

TljCp^.N)  = | 6 . (HI.38) 


It  is  to  be  emphasized  that  there  are  nine  strain  generalized 

A 

Gruneisen  parameters  for  each  mode  (p?,N)  corresponding  to  the  nine 
tensor  indices  (r,s).  For  the  isotropic  solid  the  other  tensor  elements 
may  be  expressed  for  all  longitudinal  modes  (p  ,N)  as 


li 


2C 


(11 1.39) 


and 


Y..(P0,N)  =0  for  i t j . (111.40) 

1 J X, 

Hence,  the  other  tensor  elements  of  the  strain  generalized  Gruneisen 
parameter  do  not  give  such  a simple  relation  to  the  nonlinearity  parameter 
f.  The  transverse  modes  (pt,N)  yield  a different  set  of  nine  strain 
generalized  Gruneisen  parameters. 

Using  the  expressions  for  the  longitudinal  mode  strain  generalized 
Gruneisen  parameters  [Eqs.  (II 1.37),  (1 1 1.39),  and  (1 1 1.40)],  we  obtain 
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from  r.q.  (111,13)  the  vojume  generalized  Gniueison  pnr.imct  ers  for  the 
longitudinal  modes  (pp,N)  in  terms  of  tic  clastic  constants. 


Y(pt,N)  = 


1 e - — — - 

2 3C 


(111.41) 


Using  the  relationships  among  the  second-  and  third-order  elastic 
constants  for  an  isotropic  solid  as  given  by  Toupin  and  Bernstein  (1961) 
we  get 


, 2 , ('44  + 2C144 

Y(P„.fO  = ? ! + r 


'45<>  .4  , 
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(Ill .42) 


Thus,  the  volume  generalized  Gruneisen  parameters  in  the  longitudinal 
modes  (p?,N)  contain  elastic  constants  that  are  not  measured  by  the 
harmonic  generation  technique.  However,  one  component,  ■—  B,  is  measured 
directly. 

A 

The  transverse  mode  (pt>N)  volume  generalized  Gruneisen  parameters 
arc  obtained  from  the  transverse  mode  stra i n generalized  Gruneisen 
parameters  in  the  manner  described  above  for  the  longitudinal  mode 
parameters.  They  contain  terms  involving  elastic  constants  other  than 
and  also.  Thus,  the  Gruneisen  parameter  y will  also  contain 

such  terms. 


For  an  isotropic  solid  Hq.  (TIT. 22)  simplifies  to 


YlPj.ft)  C4  + (pt  ,N)  Ct 

' Ct  + 2 Ct 


(I  II. 43) 


where  y(p^,N)  is  the  longitudinal  mode  volume  generalized  Gruneisen 
parameter  given  by  F.q.  (111.42),  y(Pt>N)  is  the  transverse  mode  volume 


generalized  Grime i sen  parameter,  and  C ^ and  C^.  are  the  weighting 
functions  given  by  liq.  (111,23). 

For  anisotropic  solids  the  calculation  of  t lie  strain  and  volume 
generalized  Gruneisen  parameters  as  well  as  the  Gruneisen  parameter  y 
is  much  more  complicated.  For  anisotropic  solids  each  longitudinal 

A 

mode  specified  by  a particular  value  of  the  set  (p^,N)  will,  in  general, 
have  nine  different  strain  generalized  Gruneisen  parameters.  This  is  in 
contrast  to  the  isotropic  solid  in  which  all  values  of  the  set  (p, ,N) 
have  the  same  set  of  nine  strain  generalized  Gruneisen  parameters.  The 
same  considerations  apply  to  the  difference  between  isotropic  ; 
anisotropic  solids  Tor  the  transverse  modes  (p  ,N). 

The  volume  generalized  Grimcisen  parameters  are  related  to  the 
st ra i n genera  1 i zed  Gruneisen  parameters  through  liq.  (111.13).  It  is 
seen  that  in  an  anisotropic  solid  the  longitudinal  and  transverse  com- 
ponents of  the  volume  generalized  Gruneisen  parameter  cannot  be 
represented  by  values  which  are  independent  of  particular  values  of 

A A 

mode  specification  (p0,N)  or  (p  ,N),  as  was  done  for  the  case  of  isotropic 
solids.  Hence,  calculation  of  the  Gruneisen  parameter  y cannot  he  done 
with  a simple  expression  as  Fq.  (111,43).  For  anisotropic  solids 
liq.  (111.22)  must  be  used.  In  addition,  for  anisotropic  solids  the 
calculation  of  the  modal  specific  heat  "weighting"  functions  |l*q.  (1 1 1. 23)] 
involves  calculation  of  directional  dependent  modal  character i s v ic  tem- 
peratures. This  further  complicates  the  problem  for  calculating 


anisotropic  Griineiscn  parameters. 


Finally,  it  is  of  interest  to  note  that  for  the  |ly0,0|  direction 
in  a cubic  crystal  the  longitudinal  mode  s c rain  generalized  Brunei  sen 
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parameter  specified  by  tensor  component  [r=l,  s=l]  may  be  written 


Yn(p£;  1,0,0)  = - 


3C  + C 

11  HI  _ 3 p 

— 2^ “ ~ p 

ZL11 


(111-44) 


This  is  identical  to  the  value  given  in  Eq.  (II  1.37)  for  al  1 
longitudinal  modes  of  an  isotropic  solid  for  the  | r= 1 , s=l]  component. 
However,  as  pointed  out  above,  Eq.  (111.44)  gives  but  one  of  an  infinite 
number  (one  for  each  direction)  of  [r=l,  s=l | strain  generalized 
Ciriineisen  parameters  which  ultimately  must  be  used  (together  with  other 
(r,s)  components)  in  the  calculation  of  the  Gruneisen  parameter  y. 


VIII.  THE  TEMPERATURE  DEPENDENCE  OF  THE 
GRUNEISEN  parameters 


It  was  assumed  in  the  expression  relating  the  thermal  expansivity 
tensor  r*  to  the  Helmholtz  free  energy  F [Eq.  (II  1.7)]  that  if  4>  (the 

J K O 

static  potential  energy  of  the  lattice  system)  and  w(p,q)  (the  frequencies 
of  the  simple  harmonic  oscillators)  were  directly  independent  of  tem- 
perature, then  the  could  be  expressed  in  terms  of  tensor  entities 
called  the  strain  generalized  Gruneisen  parameters.  These  entities  are 
also  directly  independent  of  temperature  from  their  definition 
(Eq.  (III. 9)],  and  depend  on  temperature  only  through  the  lattice 

dimensions  (i.e.,  the  y.,  arc  strain  dependent).  The  introduction  of 

J K 

the  Debye  model  eliminated  the  frequency  dependence  of  the  s t ra i n 
generalized  Griineisen  parameters  and  that  allowed  for  the  expression  of 
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the  strain  generalized  Gruneisen  parameters  exclusively  in  terms  of  the 
elastic  constants.  Hence,  the  strain  generalized  Gruneisen  parameters 
are  independent  of  temperature  to  the  extent  that  the  Debye  model  is 
accurate  and  that  the  combinations  of  elastic  constants  representing 
a particular  mode  (p,N)  and  tensor  component  (r,s)  is  independent  of 
temperature. 

It  has  been  shown  that  the  y^(p^,N)  strain  generalized  Gruneisen 
parameter  of  an  isotropic  solid  is  equal  to  y 8,  where  8 is  the  ultra- 
sonic nonlinearity  parameter.  Since  8 is  measured  directly  in  the 

A 

harmonic  generation  technique,  the  variation  of  y^j(p0,N)  with  tem- 
perature can  be  measured  directly.  Measurements  of  the  other  y (p  ,N) 

I*  S L 

A 

and  all  Yrs(Pj.,N')  for  isotropic  solids  must  be  obtained  from  other 
measurements  of  the  temperature  dependence  of  the  second-  and  third-order 
elastic  constants. 

The  volume  generalized  Gruneisen  parameters  arc  also  directly 
independent  of  temperature  since  they  involve  sums  of  the  strain 
generalized  Gruneisen  parameters  weighted  by  isothermal  compressibility 
tensor  components  No  direct  measurement  of  any  vo  1 ume  generalized 

Gruneisen  parameter  is  possible  with  the  harmonic  generation  technique 
since  all  such  parameters  involve  terms  other  than  C and  Gj ^ (isotropic 
solids). 

Since  the  volume  generalized  Gruneisen  parameters  arc  directly 


independent  of  temperature,  the  Gruneisen  parameter  y can  depend 
directly  on  temperature  only  through  the  modal  specific  heats  ns  shown 
in  Eq.  (111.22)  [Eq . (I  [1.43)  for  isotropic  solids].  The  modal  specific 
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heats  can  be  calculated  (usually  numerically)  as  functions  of  tem- 
perature and  the  temperature  dependence  of  all  the  elastic  constants  can 
be  measured  experimentally  by  combining  the  harmonic  generation 
technique  with  other  techniques  [for  example,  an  ultrasonic  beam  mixing 
technique  (Dunham  and  Huntington,  1970)].  Hence,  all  the  factor s in  the 
right  hand  side  cf  F.qs.  (111.22)  or  (111,43)  can  be  measured  or  calcu- 
lated as  functions  of  temperature,  and  the  Crtlneiscn  parameter 
calculated. 

The  Grli.. eisen  parameter  can  also  be  determined  experimentally  from 
other  measurements  [Eq.  (111,1)] , Comparing  the  theoretically  determined 
value  for  the  Gruneisen  parameter  [Eqs.  (IIT.22)  or  (1 1 1.43)  and  the 
experimentally  determined  value  [Hq.  (TI1.1)]  gives  a direct  check  of 
the  validity  of  the  Debye  theory  for  a particular  solid. 

Although  some  comparisons  between  the  theoretically  and  experimentally 
determined  Crtlneiscn  parameters  have  been  made  for  various  solids  (Cerlich, 
1975;  Bruggcr  and  Fritz,  1967;  Schucle  and  .Smith,  1963;  Sheard,  1958),  the 
results  are  largely  inconclusive  because  of  the  lack  of  experimental  data 
on  the  temperature  dependence  of  the  third-order  elastic  constants.  This 
lack  of  information  provided  motivation  for  the  measurements  of  the 
temperature  dependence  of  third-order  elastic  constants  reported  in  this 
d issertat ion . 


CHAPTER  IV 
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APPARATUS,  PROCEDURES,  AND  SAMPLES 
I.  EXPERIMENTAL  CONSIDERATIONS 

The  determination  of  the  adiabatic  nonlinearity  parameters  of  the 
solid  samples  was  made  by  measuring  both  the  amplitudes  of  the 
longitudinal  wave  introduced  into  the  specimen  and  the  amplitude  of  the 
generated  second  harmonic  component.  The  fundamental  frequency  was 
chosen  to  he  30  MHz.  Such  a choice  represents  a compromise  between  the 
fact  that  the  second  harmonic  amplitude  A.,  increases  as  the  square  of 
the  frequency  and  the  fact  that  attenuation  and  effects  of  nonpara  1 lei i sm 
of  the  sample  surfaces  also  increases  with  frequency.  To  avoid 
diffraction  effects  it  was  desirable  to  work  above  10  Mill  (Gaustcr, 

1966).  for  the  samples  used  in  these  measurements  a »0  MHz  fundamental 
frequency  provided  an  adequate  signal -to-noi se  ratio  while  introducing 
errors  due  to  attenuation  that  were  negligible  relative  to  other  sources 
of  error. 

A pulse-echo  technique  was  used  in  all  measurements  to  eliminate 
complicated  interference  effects.  A 60  Ilz  pulse  repetition  rate  was 
chosen  to  minimize  heating  of  the  samples  and  to  reduce  ground  loops. 

Since  only  in  the  limit  of  infinitesimal  fundamental  amplitude  is 
the  second  harmonic  related  in  amplitude  to  the  fundamental  according 
to  a simple  integral  power  law  (Thurston  and  Shapiro,  1967),  the  funda- 
mental amplitude  was  kept  as  low  as  possible  consistent  with  a useful 
signal -to-noi sc  ratio. 


! 


I 


32 


II.  ABSOLUTE  AMPLITUDE  MEASUREMENTS  OF  Till' 
ULTRASONIC  WAVE  HARMONICS 
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The  amplitudes  of  the  fundamental  and  generated  second  harmonic 
waves  were  measured  absolutely  in  each  sample  at  room  temperature.  At 
all  other  temperature  points  the  amplitudes  were  measured  relative  to 
the  room  temperature  values.  Irom  a knowledge  of  the  absolute  wave 
amplitudes,  the  third-order  elastic  constant  C^  is  calculated  from 
Eq.  (11.23)  as 


C 
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(IV.  1) 


where  ^ is  calculated  from  velocity  measurements. 

Room  Temperature  Apparatus 

The  absolute  amplitude  measurements  were  made  with  the  room 
temperature  apparatus  as  shown  in  Figure  IV. 1.  A cross-sect  tonal  view 
of  the  apparatus  is  shown  in  Figure  IV. 2.  An  ultrasonic  signal  is 
generated  by  an  X-cut  piezoelectric  quartz  transducer,  propagates 
through  the  sample,  and  is  detected  at  the  other  end  by  an  air-gap 
capacitive  receiver.  The  capacitive  receiver  is  a parallel  plate 
capacitor  consisting  of  an  electrically  conducting  sample  face  and  0.400 
inch  diameter  electrode  recessed  approximately  7 microns  from  the  sample 
surface.  A dc  bias  of  the  order  of  140  volts  is  applied  across  the  gap 
through  a large  resistor  (approximately  1 MO)  to  keep  the  current 
through  the  detector  low  in  case  of  electrical  arc  over. 

Tn  order  to  make  accurate  amplitude  measurements  of  the  ultrasonic 
waves,  provision  is  made  to  place  a substitutional  signal  on  the  receiver 
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through  a large  resistor  (approximately  10  Kft)  between  the  receiver  and 
the  substitutional  signal  source. 

Room  Temperature  Calibration  and  Procedure 

The  absolute  determination  of  the  fundamental  and  second  harmonic 
ultrasonic  wave  amplitudes  is  made  with  the  apparatus  arrangement  shown 
in  the  block  diagram  of  Figure  IV. 3. 

A cw  signal  from  a stable  VFO  is  used  to  drive  a gated  amplifier. 

The  resulting  pulsed  30  MUz  signal  is  passed  through  a power  amplifier 
designed  to  deliver  full  power  at  any  load  impedance.  The  output  from 
the  power  amplifier  is  then  passed  through  a 30  MHz  bandpass  filter 
which  insures  the  spectral  purity  of  the  ultrasonic  wave,  even  when  the 
transducer  is  driven  off  resonance. 

At  the  opposite  end  of  the  sample  the  ultrasonic  wave  generates  an 
electrical  signal  in  the  capacitive  receiver.  This  signal  is  sent 
through  cither  a 30  MHz  or  60  MHz  bandpass  amplifier.  From  there  the 
signal  is  taken  to  a boxcar  integrator.  The  boxcar  integrator  selects 
a portion  of  the  first  echo  and  gives  an  output  signal  proportional  to 
the  time  average  of  that  portion  of  the  echo. 

After  recording  the  magnitude  of  the  fundamental  and  second  harmonic 
outputs, a continuous  (cw)  substitutional  signal  is  introduced  at  the 
capacitive  receiver.  The  30  MHz  filter  used  previously  to  insure  the 
spectral  purity  of  the  gated  signal  is  now  used  to  filter  the  30  MHz  cw 
signal.  The  60  MHz  cw  substitutional  signal  is  derived  by  passing  the 
30  MHz  cw  signal  through  a ring  bridge  mixer  and  filtering  the  output 
with  a 60  MHz  bandpass  filter.  The  magnitudes  of  the  two  signals  arc 
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igure  IV. 3.  Block  Diagram  for  the  Room  Temperature  Nonlinearity  Measurements 
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adjusted  to  give  the  same  output  signals  as  recorded  in  the  pulse 
measurements.  The  cw  signals  arc  then  measured  with  an  rf  voltmeter, 
and  the  amplitudes  of  the  fundamental  and  second  harmonic  wave  components 
are  calculated. 

The  equivalent  circuit  for  making  the  ultrasonic  wave  amplitude 
calculations  is  shown  in  Figure  JV.4.  A Norton  equivalent  circuit  is 
used  to  represent  the  ac  equivalent  capacitive  receiver  circuit.  The 
"switch"  Sj  is  opened  by  turning  off  the  ultrasonic  signal  to  the 
receiver.  "Switch"  S?  is  closed  by  connecting  a signal  generator  to  the 
apparatus. 

It  has  been  shown  (Peters,  1968)  that  the  Thevenin  equivalent 
circuit  of  the  capacitive  receiver  is  given  as  a voltage 

2 AV, 

V = — (IV.  2) 
b0 


feeding  the  load  resistance  R (amplifier  impedance)  through  a capacitance 
of  value  C^t  the  quiescent  receiver  capacitance.  In  the  above  equation, 

A is  the  amplitude  of  the  ultrasonic  wave  under  consideration,  is  the 
capacitive  receiver  bias  voltage  (140  volts  in  these  experiments),  and 
Sg  is  the  receiver  gap  spacing. 

Since  the  receiver  voltage  is  related  to  the  receiver  current  by 
iD  = VwCp  where  « is  the  angular  frequency  of  the  ultrasonic  wave,  the 
Norton  equivalent  circuit  (Bains,  1974)  of  the  detector  is  a current 


2AVS 


(TV. 3) 


in  parallel  with  the  resistance  R and  receiver  capacitance  as  shown 
in  Figure  IV. 4.  The  substitutional  signal  source  is  represented  in 
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Figure  IV. 4 as  the  current  generator  Gg.  When  the  magnitude  of  this 
current  source  is  adjusted  to  give  the  same  output  from  the  amplifier  as 
the  pulsed  ultrasonic  signal,  ig  is  equal  to  i.^.  It  is  then  a matter  of 
measuring  ig  to  find  the  ultra  sonic  wave  amplitude  from  Eq.  (TV.  S'). 

Since  voltages  are  generally  much  easier  to  measure  than  currents, 
ig  is  calculated  from  measurement  of  the  substitutional  source  voltage 
and  knowledge  of  the  impedance  through  which  the  current  flows.  The 
inductance  L shown  in  Figure  IV. 4 is  the  inductance  of  the  wire  leading 
from  the  detector  to  the  BNC  connector  (see  Figure  IV. 2,  page  35)  and 
the  capacitance  Cg  is  the  lumped  stray  capacitance  of  the  receiver. 

At  the  frequencies  used  in  these  experiments  the  resistor  R between 
the  substitutional  source  and  the  receiver  does  not  behave  as  a pure 
resistance.  Hence,  the  impedance  of  the  resistor  must  he  measured  at 
each  frequency  used.  The  equivalent  circuit  for  measurement  of  the 
resistor  impedance  is  shown  in  Figure  IV. 5. 

The  sample,  bottom  plate,  and  detector  assembly  are  removed  from 
the  apparatus  and  50 ft  terminators  are  placed  on  the  ’’signal  in”  and 
’’signal  out”  BNC  connectors.  A cw  signal  generator  is  connected  to  the 
"signal  in”  terminator  designated  as  point  1 in  Figure  IV. 5.  Probes  A 
and  B (with  isolator  tips)  of  a vector  voltmeter  (Hewlett-Packard 
Model  8405A)  are  placed  at  point  1,  the  phase  angle  between  the  signals 
entering  the  probes  are  zeroed  and  the  signal  amplitudes  arc  measured. 
Probe  B of  the  vector  voltmeter  is  then  placed  at  point  2 while  leaving 
probe  A at  point  1.  The  signal  generator  is  readjusted  to  give  the 
original  reading  for  channel  A and  the  amplitude  of  the  B channel  and 
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the  phase  angle  between  the  channels  are  recorded.  The  impedance  of  the 
resistor  R is  then  given  by 


[ jwC  + 


-1  VB1 


B2 


jwL 


IV. 4) 


where  VD1  is  the  signal  amplitude  measured  by  probe  B at  point  1,  V’  is 

the  amplitude  measured  by  probe  B at  point  2,  and  <p  is  the  phase  angle 

betweei.  Vn1  and  V . L is  the  inductance  mentioned  previously,  w is  the 

angular  frequency  of  the  cw  input  signal,  Rj  is  the  measured  resistance 

of  the  "signal  out"  50ft  terminator,  C is  the  stray  capacitance  at  point  2 

2 

(including  the  probe  tip  capacitance),  and  j = -1. 

The  substitutional  current  is  calculated  from 


s " 2 + + V + jo-rj^rl 

where  is  the  voltage  across  the  current  generator  (Figure  IV. 4, 
page  59) . 

III.  RELATIVE  AMPLITUDE  MEASUREMENTS  OF  THE 
ULTRASONIC  WAVE  HARMONICS 


The  Cryogenic  Apparatus 

The  apparatus  used  to  make  the  relative  ultrasonic  wave  amplitude 
measurements  basically  consists  of  a transducer  driver  assembly,  a sample 
holder,  and  an  air-gap  capacitive  receiver  enclosed  in  a stainless  steel 
can  as  shown  in  Figure  IV. 0.  This  can  is  surrounded  by  another  can  so 


that  the  space  between  the  cans  can  he  evacuated.  This  provides  an 
insulating  jacket  around  the  inner  can.  The  cans  are  supported  by  three 
thin-walled  cupro-nickcl  tubes.  Two  of  the  tubes  have  a smaller  cupro- 
nickel tube  centered  inside  the  diameters  of  which  were  chosen  to  provide 
a 50ft  coaxial  transmission  line:  one  for  rf  input  to  the  transducer  and 

the  other  for  the  output  from  the  detector.  These  tubes  are  also  used 
as  vacuum  lines  to  the  detector  cavity  and  the  outer  can.  The  other  tube 
houses  the  leads  to  the  temperature  sensors  and  to  the  heater,  and  also 
is  used  to  supply  a pressure  to  the  inner  can. 

The  gap  spacing  is  controlled  pneumatical ly  as  follows.  The 
capacitive  receiver  is  similar  to  that  of  the  room  temperature  apparatus, 
except  that  the  lapped  ring  against  which  the  sample  rests  is  undercut 
approximately  .015  inch  to  make  it  a flexible  diaphragm  as  shown  in 
Figure  IV. 7.  The  space  above  the  sample  is  evacuated  (through  the  tube 
which  is  also  the  transmission  line  for  the  signal  from  the  detector) 
and  the  rest  of  the  sample. 

The  entire  apparatus  is  placed  inside  a helium  research  dewar. 

After  sealing  the  dewar  and  pumping  on  the  coolant  (either  liquid 
nitrogen  or  liquid  helium),  temperature  below  the  room-pressure  boiling 
point  of  the  coolant  may  be  obtained. 

The  temperature  around  the  sample  is  controlled  by  an  electrical 
resistance  heater  which  is  connected  to  a standard  commercial  temperature 
controller. 

The  temperature  measurements  are  made  by  sending  a measured  current 
from  a constant  current  generator  through  the  sensors  (platinum  or 
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germanium)  mounted  next  to  the  sample  and  measuring  the  voltage  across 
the  sensors  with  a high  input  impedance  (1000  MO)  dc  voltmeter.  The 
resistance  of  the  sensor  is  calculated  and  the  temperature  is  determined 
from  a sensor  resistance  versus  temperature  calibration  curve  furnished 
by  the  manufacturer. 


C ryogen i c Non  1 i nca r i tv  Measurements 


The  relative  amplitude  cryogenic  measurements  are  made  with  the 
equipment  arrangement  shown  in  figure  IV. 8.  The  signal  paths  are  the 
same  as  with  the  room  temperature  apparatus  except  that  no  substitutional 
signal  is  introduced  here. 

The  basis  for  the  experimental  procedure  for  making  relative 

amplitude  measurements  may  he  seen  from  an  examination  of  the  expression 

for  calculating  the  third-order  clastic  constant  fjjj  Uq*  IV.  1) 

8 ^7 

| C 1 1 j = - — y -7- ■ + 1)3  Cjj].  It  is  seen  that  if  the  fundamental  wave 

A “k~a 

amplitude  Aj  is  maintained  constant  as  a function  of  temperature,  then 
changes  in  result  from  changes  in  Cjj  and  the  second  harmonic 

amplitude  A.,.  Changes  in  (! ^ are  determined  from  independent  measurements, 
To  maintain  Aj  constant,  it  is  necessary  only  to  maintain  a constant 
30  MHz  amplifier  output  into  the  boxcar  integrator  by  suitable  changes  in 
the  rf  drive  voltage  the  transducer.  This  is  done  while  maintaining 
a constant  bias  voltage  across  the  capacitive  receiver  and  keeping  the 
receiver  spacing  (and  hence  capacitance)  constant. 

With  Aj  held  constant  as  a function  of  temperature,  the  relative 
measurements  are  made  on  the  second  harmonic  amplitudes.  Let  the  output 
of  the  60  MHz  amplifier  be  expressed  as  V . Let  the  amplitude  of  the 


4K 


second  harmonic  be  A7  and  the  capacitive  receiver  bias  voltage  I : V^, 
If  the  capacitance  of  the  capacitive  receiver  is  held  constant  as  a 
function  of  temperature,  then  so  is  the  gap  spacing.  Hence,  we  may 
write  from  fiq.  (IV. 2)  that  at  some  temperature  T 


V0  = KA2Vb 


(IV. h) 


where  K is  a constant.  At  the  reference  temperature  T (room 
temperature),  we  may  write  the  equation  as 


V R - KA  RV  R 
V0  - KA  Vb  . 


(IV. 7) 


If  the  bias  voltage  is  altered  at  temperature  T so  that  the  output 
of  the  60  MHz  amplifier  remains  constant,  i.c.. 


V = V 
0 0 


(IV. 8) 
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Hence,  the  relative  amplitude  measurements  of  the  second  harmonic  may  be 
obtained  from  measurements  of  t he  capacitive  receiver  bias  voltage. 

It  can  be  shown  from  expression  (IV.9)  and  Eq.  (IV. 1)  that  the 
ratio  of  at  temperature  T to  CR^  at  temperature  TR  is  (Peters, 


R V R C 
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(IV. 10) 
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where 


flR  - * 
6 - 3 


A 212 

A j k a 


(IV. 11) 


k = 1 + o(T  - TK) 


(IV. 12) 


and  a is  the  thermal  expansion  coefficient,  which  for  fused  silica  is 
less  than  10  (> , which  is  negligible  in  these  experiments. 


IV.  VELOCITY  MEASURr.MF.NTS 

The  second-order  clastic  (SOI:)  constants  Cjj  appear  in  both  Eqs. 

(IV.  1)  and  (IV.  10).  is  calculated  from  the  expression 

Cu  = pV2  (IV. 13) 

whore  p is  the  mass  density  of  the  material  and  V is  the  ultrasonic  wave 
velocity  in  the  appropriate  direction.  For  Raised  silica  the  sonic 
velocity  is  independent  of  direction.  The  SOI:  constants  were  determined 
absolutely  from  room  temperature  velocity  measurements  by  an  ultrasonic 
pulse-echo  technique.  The  tabulated  data  of  A.  Zarcmbowitch  (1976)  were 
used  in  the  calculation  of  the  relative  cryogenic  SOF.  constants. 

The  Capacitive  Driver-Capacitive  Receiver  System 

Until  recently,  it  has  been  necessary  to  generate  the  ultrasonic 
pulses  and  to  detect  the  resulting  echoes  in  velocity  measurements  by 
applying  bursts  of  rf  to  a piezoelectric  transducer  bonded  to  the  specimen 
surface.  Such  an  arrangement  necessitates  making  a tedious  correction 
in  the  velocity  measurement  due  to  ultrasonic  wave  phase  shifts  at  the 


so 


specimen -bonded  transducer  interface.  In  sone  situations,  the  effect  of 
this  coupling  between  specimen  and  transducer  leads  to  considerable 
error  if  neglected  (McSkimin  and  Andreatch,  1962). 

In  order  to  eliminate  the  bond  problem  in  compressional  wave 
measurements,  the  air-gap  capacitive-receiver  used  previously  to  measure 
amplitudes  of  ultrasonic  waves  (Gauster  and  Breazeale,  1966;  Peters, 
Breazeale,  and  Par6,  1968)  was  modified  to  function  as  an  air-gap 
capacitive  driver  as  well  (Cantrell  and  Breazeale,  1974).  The  air-gap 
capacitive  driver-capacitive  receiver  combination  allowed  the  sample  ends 
to  vibrate  with  free-free  boundary  conditions  exactly.  No  bond 
corrections  were  necessary  because  there  were  no  bonds. 

More  recently,  a dielectric  capacitive  driver  was  developed  to 
replace  the  air-gap  driver  (Cantrell  and  Breazeale,  1977).  The  new  driver 
is  capable  of  generat  ing  much  larger  longitudinal  ultrasonic  pul  sos  than  the 
previous  one  because  one  can  apply  larger  voltages.  I;rec-free  boundary 
conditions  stil  1 exist  at  t lie  sample  surfaces  with  this  driver  and  the 
error  in  velocity  measurements,  without  making  any  corrections,  is  com- 
parable to  the  error  in  velocity  measurements  made  with  bonded  transducers, 
after  corrections  for  bond  phase  shifts. 

The  mechanical  parts  of  the  dielectric  capacitive  driver-capaci tive- 
receiver  system  are  the  same  as  the  room  temperature  apparatus  (see 
figure  IV. 2,  page  35),  except  that  the  substitutional  signal  capabilities 
are  not  used  and  the  quartz  transducer  is  replaced  with  a piece  of  mica 
0.625  inch  in  diameter  and  approximately  0.0002  inch  thick.  The  mica  is 
bonded  to  neither  the  sample  nor  the  driver  button. 


A schematic  of  the  dielectric  capacitive  driver  is  shown  in 
Figure  IV. 9.  The  sample  surface  (outer  electrode)  is  designated  as 
(1)  and  the  driver  button  (inner  electrode)  is  designated  as  (5). 

Between  the  electrodes  and  the  central  mica  dielectric  (5),  two  air 
layers  (2  and  4)  are  present  because  of  slight  irregularities  of  the 
contact  surfaces. 

A pulsed  sinusoidally  varying  rf  voltage,  V = V^sinmt,  applied 
between  the  electrodes  (1  and  5),  produces  sinusoidally  varying  electric 
fields  I.j,  Ii9,  and  lu.  It  has  been  demonstrated  (I.egros,  Lcwincr,  and 
Biquard,  1972)  that  the  generation  of  ultrasonic  waves  in  this  situation 
is  due  to  the  effect  of  electrodynamic  forces  acting  directly  on  the 
ultrasonic  wave  propagation  medium.  Hence,  if  F.^  is  the  electric  field 
near  the  surface  of  the  sample,  the  force  per  unit  area,  p,  at  this 
surface  is 


1 r 2 
P = 2 eoE3 


where  is  the  dielectric  permittivity  of  air.  Since 


(IV. 14) 
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sintot 
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(IV. 15) 


1 e E_  ^sin^o>t  = ~ z E ^ (1  - cos  2wt) 

2 o 3o  4 o 3o 


(IV. 16) 


For  present  purposes,  it  is  adequate  to  ignore  the  constant  term  in 
Eq.  (IV. 16)  and  consider  only  the  cos  2wt  term.  One  sees  that  the 
pressure  on  the  sample  surface  is  applied  at  twice  the  frequency  of  the 
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rf  voltage*  Thus,  the  ultrasonic  wave  generated  in  the  sample  has  twice 
the  frequency  of  the  applied  rf  voltage. 

A Tuned  Transmission  Line  Voltage  Transformer 

The  capacitive  transducer  is  driven  with  a pulsed  rf  signal.  The 
frequency  is  varied  around  15  MHz  to  produce  successive  resonances.  To 
obtain  the  voltages  necessary  for  adequate  ultrasonic  amplitudes,  a 
specific  length  of  transmission  line  was  used  to  *■  as  a tuned  voltage 
transformer  between  the  gated  amplifier  and  the  capacitive  driver. 

If  S is  the  length  of  the  transmission  line,  f the  rf  drive 
frequency,  the  voltage  across  a purely  reactive  load  of  magnitude 
]Zj|  (capacitive  driver),  the  voltage  out  of  the  gated  amplifier, 
and  the  characteristic  impedance  of  the  transmission  line,  then  for 
a lossless  transmission  line 

Vc  = Vn  sin*  esc (*  - — ) (IV. 17) 

o U C 

where 

-i  lzJ 

<J>  = - tan  - (IV.  18) 

i) 

w = 2irf  (IV. 19) 

and  c is  the  velocity  of  propagation  of  electromagnetic  waves  in  the 

VS 

transmission  line.  Thus,  for  the  proper  length  S,  — becomes  very  large. 

V0 

But  the  impedance  looking  out  of  the  gated  amplifier  into  the  trans- 
mission line  behaves  as 

Zin  = -jZ„  ta"(*  - f3  • 


(IV. 20) 
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Thus,  for  the  condition  that  ^—becomes  large,  Z.  , and  consequently, 

0 

V^,  become  small.  For  this  reason,  the  length  of  transmission  line  is 
chosen  such  that  looks  inductive  for  our  operation  frequency.  A 
variable  capacitor  is  then  placed  in  parallel  with  to  form  a 
parallel  resonant  circuit.  This  arrangement  increases  the  intensity  of 
the  generated  ultrasonic  signal  by  approximately  30  dB  over  that  of  the 
untuned  transmission  line. 

The  primary  limitation  on  the  generated  ultrasonic  wave  amplitude 

is  set  by  dielectric  breakdown  in  the  capacitive  transducer.  However, 

this  limitation  is  not  as  serious  as  one  might  assume.  As  the  spacing 

of  the  electrodes  in  air  is  decreased,  the  breakdown  potential  goes  up 

almost  exponentially  (Cantrell  and  Breazeale,  1974).  In  the  air-gap  driver  we 

KV 

were  able  to  use  fields  as  large  as  200  — without  breakdown.  With  the  mica 

v cm 

dielectric  driver,  we  nave  been  able  to  use  fields  larger  than  1000 
Method  of  Measurement 

The  velocity  measurements  were  made  with  the  equipment  arrangement 
as  shown  in  Figure  IV. 10.  The  gated  double-pulse  superposition  method 
of  Williams  and  Lamb  (1958)  was  used.  A pulsed  ultrasonic  wave  is  fed 
into  the  sample  from  the  dielectric  capacitive  driver  followed  by  another 
pulsed  ultrasonic  signal  phase-locked  but  delayed  in  time  with  respect 
I to  the  first.  The  delay  is  adjusted  such  that  superposition  of  the 

desired  echoes  from  the  two  pulse  trains  is  achieved.  The  resulting 
signal  is  received  at  the  other  end  of  the  sample  by  the  air-gap 
' capacitive  receiver,  passes  through  a 30  MHz  bandpass  amplifier  and  is 

displayed  on  an  oscilloscope.  Utilizing  the  fact  that  a continuous 

L. 
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Figure  IV. 10.  Block  Diagram  of  Experimental  Apparatus  for  Velocity 
Measurements . 
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change  in  ultrasonic  frequency  results  in  a succession  of  superimposed 
echo  maxima  and  minima,  the  ultrasonic  velocity  is  calculated  from 

v = 2nP  ^ (IV. 21) 

An 

where  Af  is  the  change  in  ultrasonic  frequency  corresponding  to  An 
minima,  l is  the  sample  length,  and  m is  an  integer  (2  in  these  experi- 
ments) giving  the  relative  position  of  the  initial  pulse  train  and  the 
delayed  pulse  train.  To  enhance  the  sensitivity  in  detecting  the  super- 
imposed echo  minima,  a boxcar  integrator  was  incorporated  into  the  system 
between  the  30  MHz  amplifier  output  and  the  oscilloscope. 


V.  SAMPLES  AND  SAMPLE  PREPARATION 

Samples 

The  fused  silica  samples  used  in  this  investigation  were  obtained 
from  Amersil,  Inc.  of  Hillside,  New  Jersey,  and  were  manufactured  under 
the  commercial  designation  Suprasil.  Suprasil  represents  a particular 
grade  of  fused  silica  which  has  been  given  a rating  of  Grade  A, 

Type  458-677,  according  to  military  specification  MIL-G-174. 

Sample  homogeneity  is  measured  by  the  maximum  variation  in  mass 
density  of  the  material  and  in  optics  is  most  generally  expressed  as  a 
change  in  refractive  index  (An)  along  an  optical  path.  The  Suprasil 
samples  arc  designated  as  optical  grades  1 and  2.  The  designation  1 
implies  a stated  maximum  variation  of  refractive  index  in  all  spatial 
directions,  whereas  the  designation  2 implies  a stated  maximum  refractive 
index  variation  in  only  one  direction  (perpendicular  to  the  sample 


surfaces) . 
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Tlie  designation  W is  used  to  distinguish  the  OH  content  of  the 
samples.  Those  samples  represented  by  the  designation  W contain  approxi- 
mately 5 ppm  of  OH  content,  whereas  the  unlettered  samples  have  OH 
contents  of  approximately  1200  ppm. 

Although  all  Suprasil  samples  contain  less  than  1 ppm  of  total 
metallic  impurities,  the  Suprasil  W samples  contain  approximately  260  ppm 
each  of  fluorine  and  chlorine  impurities,  whereas  the  unlettered  samples 
contain  about  130  ppm  each  of  fluorine  and  chlorine. 

All  samples  used  in  these  experiments  were  approximately  one  inch 
in  diameter.  The  length  of  each  sample,  the  above  mentioned  impurity 
properties,  and  the  measured  ultrasonic  velocities,  at  room  temperature 
(27  °C),  are  tabulated  in  Table  IV. 1. 


TABLE  IV. 1 

SAMPLES  USED  IN  THIS  STUDY  AND  THEIR  PROPERTIES 


Sample 

Impurity  Content  (ppm) 
OH  Fluorine  Chlorine 

Sample  Length 
(cm) 

Longitudinal  Ultrasonic 
Velocity  (xlO^  cm/sec) 

Suprasil — W 1 

5 

260 

260 

1 . 2590±0. 0003 

5.952610.0016 

Suprasil — W2 

5 

260 

260 

1.2587+0.0002 

5.953310.0011 

Suprasi 1 1 

1200 

130 

130 

1 . 2593±0. 0002 

5.936310.0012 

Suprasil  2 

1200 

130 

130 

1.2586+0.0002 

5.936310.0011 

Since  the  fused  silica  samples  are  nonconductive,  a conductive  copper 

o 

surface  of  the  order  of  1000  A thick  was  evaporated  onto  the  ends  of  the 
samples  to  form  the  ground  potential  electrodes  for  the  transducer  and 


detector. 


The  ultrasonic  velocity  measurements  of  the  samples  were  made  with 
the  capacitive  driver-capacitive  receiver  system  and,  hence,  no  trans- 
ducers were  bonded  to  the  sample  surfaces  for  these  measurements.  The 
use  of  a capacitive  driver-capacitive  receiver  system  for  the  low- 
temperature  relative  measurements  of  C^  would  require  extensive  modifi- 
cation of  the  cryogenic  apparatus  to  produce  fundamental  ultrasonic  wave 
amplitudes  large  enough  to  generate  measurable  second  harmonics.  Further, 
the  amplitude  measurements  of  the  first  and  second  harmonics  do  not 
depend  upon  superposition  of  ultrasonic  pulse-ccho  trains  as  did  the 
velocity  measurements.  Hence,  corrections  for  transducer  bonds  arc  not 
important  for  these  measurements.  For  expediency,  then,  piezoelectric 
transducers  wore  bonded  to  the  sample  surfaces  for  the  nonlinearity 
measurements . 

The  transducers  used  in  the  nonlinearity  measurements  were  X-cut 
quartz  crystals  bonded  to  the  sample  with  the  adhesive  from  a Commercially 
available  cellulose  tape  (manufactured  by  the  Technical  Tape  Corporation, 
Carbondale,  Illinois,  under  the  designation  7510-551-9818,  L590C,  Type  1, 
Class  A) . 

The  adhesive  was  applied  in  the  manner  described  by  Yost  (1972).  A 
hole  was  cut  in  a plastic  film  having  a diameter  slightly  larger  than  the 
diameter  of  the  transducer.  The  hole  was  centered  on  the  sample  surface 
and  a strip  of  tape  was  put  over  the  hole.  The  tape  stuck  to  the  sample 
through  the  hole  and  air  bubbles  were  rubbed  out  with  a stick.  IVater  was 
then  applied  to  the  tape  backing  and  left  until  the  water  turned  cloudy 
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(I  to  5 minutes).  The  water  was  then  blotted  away  and  the  tape  was 
gently  lifted  off.  An  amount  of  adhesive  covering  an  area  determined 
by  the  hole  was  left  on  the  sample  surface.  The  sample  was  left  in  a 
dessicator  for  twenty-four  hours  to  remove  excess  moisture  in  the 
adhesive.  The  quartz  transducer  was  then  applied  to  the  adhesive  sur- 
face and  bonded  by  rubbing  the  transducer  firmly  with  the  end  of  a 
teflon  cylinder. 


CHAPTER  V 


RESULTS  AND  DISCUSSION 

I.  ROOM  TEMPERATURE  NONLINEARITY  MEASUREMENTS 


The  results  of  the  absolute  amplitude  measurements  at  room 

temperature  for  all  the  fused  silica  samples  are  given  in  Tabic  V.l. 

2 

A plot  of  A0  versus  A^  for  each  sample  is  shown  in  Figures  V.l,  V.2, 
V.3,  and  V.4.  The  straight  line  in  each  of  these  figures  is  a least 
squares  fit  to  the  data  points.  It  is  seen  that  t ho  data  are  fit  very 
well  with  a straight  line.  The  lines  do  not  pass  through  the  origin  and 
the  resulting  intercept  gives  some  measure  of  residual  electronic 
noise  (Bains,  1974)  and  is  seen  to  be  small  for  these  measurements. 

The  slope  of  the  line  is  (in  the  least  squares  sense)  the  best  fit 
to  the  equation 

a2  = | s (\.n 

and  may  be  used  to  calculate  the  nonlinearity  parameter  8.  However,  a 
much  more  sensitive  method  to  calculate  8 was  developed  by  Yost  (1972). 
In  his  method  8 is  calculated  from 


*The  measurements  for  the  Suprasil  1 sample  were  repeated  after 
noticing  that  the  An/A  “ values  of  this  sample  were  unusually  high  compared 
to  the  A2/Aj-  values  of  the  other  samples.  The  repeated  data  confirms 
the  original  measurements  and  eliminates  the  possibility  that  electronic 
instability  introduced  the  high  values  of  A^/A^. 
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(V.2) 


and  is  plotted  as  a function  of  Aj.  This  method  allows  measurements 
with  good  signal -to-noise  ratios,  and  extrapolating  the  curve  to  zero 
value  of  amplitude  gives  a value  of  0 which  meets  the  theoretical 
assumption  of  infinitesimal  amplitude  waves  (Thurston  and  Shapiro,  1967). 

The  B versus  Aj  curves  for  each  sample  is  given  in  Figures  V.5,  V. 6, 
V.7,  and  V. 8 and  include  the  corrections  for  electronic  noise  obtained 
from  the  curves  of  Figures  V.l,  V.2,  V.3,  and  V.4.  The  extrapolated  0 
values  arc  listed  in  Table  V.2,  along  with  calculated  values  of  Cj j and 
Cjjl  at  room  temperature  for  each  of  the  samples  used  in  these  experiments 
The  second-order  elastic  constants  are  calculated  from 


Cn  = P0V-  CV.3) 


where  pQ  is  the  mass  density  of  the  samples  and  V the  longitudinal 
ultrasonic  velocity.  The  third-order  elastic  constants  are  calcu- 

lated from 


cm  = -3Cn(0  ♦ 1)  . 


II.  TEMPERATURE  DEPENDENCE  OF  0 AND  C 


(V,  4 ) 


111 


Table  V.  3 gives  the  measured  relative  values  of  3 and  as  a 

function  of  temperature.  The  temperature  dependence  of  Cjj  was  calcu- 
lated from  the  data  of  A.  Zarembowitch  (1976)  who  measured  as  a 
function  of  temperature  for  a sample  of  Puropsil  fused  silica.  The 
temperature  dependence  of  the  samples  used  in  this  experiment  are 
assumed  to  have  the  same  relative  temperature  dependence  as  that  of 
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TABLE  V. 3 

MEASURED  RELATIVE  VALUES  OF  3 AND  C 

AS  A FUNCTION  OF  TEMPERATURE  111 


Samp  1 c 

T (°K) 

lei 

1 2 

C 1 1 1 (10  “ dynes/ cm"") 

Suprasil  W1 

2.9 

3.765 

6.276 

5 . 4 

3.765 

6.271 

6.5 

3.763 

6.262 

8.9 

3.762 

6.254 

18.8 

3.768 

6.256 

30 

3.709 

6.107 

42 

3.746 

6.178 

53 

3.755 

6.166 

65 

3.810 

6.307 

77 

3.743 

6.  162 

88 

3.584 

5.814 

100 

3 . 696 

6.070 

112 

3.744 

6.187 

123 

3.723 

6.154 

136 

3.732 

6.184 

147 

3.823 

6.403 

159 

3.760 

ro 

cn 

170 

5.671 

6 . 090 

181 

3.815 

6.4  32 

194 

3.824 

6.470 

205 

3.842 

6.522 

217 

5.794 

6.427 

233 

3.812 

6.490 

245 

3.824 

6.535 

259 

3.825 

6.556 

273 

3.842 

6.618 

282 

co 

cn 

Cn 

6. 663 

295 

3.859 

6.695 

3.0 

3.732 

6.199 

5.1 

3.803 

6.353 

10.1 

3.975 

6.726 

18.8 

3.697 

6.093 

26.2 

3.766 

6.236 

36 

3. 846 

6.405 

45 

3.895 

6.506 

77 

3.811 

6.310 

88 

3.827 

6.353 

100 

3.785 

6.265 

112 

3.783 

6.270 

123 

3.681 

6.053 

136 

3.762 

6.247 

147 

3. 763 

6.262 

Suprasil  IV 2 


gure  V.9.  B versus  Temperature  for  Suprasil  Wl. 


Figure  V.13.  (iraph  of  C versus  Temperature  for  Supras 


water  content  in  the  bonding  material  (cellulose  tape  glue).  A 
reasonable  estimate  would  be  that  the  random  error  contributed  by 
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temperature  changes  in  the  transducer  bonds  is  no  more  than  ± 3 percent. 
This  number  was  obtained  from  a comparison  of  measurements  at  a given 
temperature  after  cycling  through  a temperature  range  of  approximately 
120  degrees. 

In  summary,  the  total  error  in  the  room  temperature  measurements 
of  the  absolute  ultrasonic  wave  amplitudes  is  estimated  to  he  no  greater 
than  + 14  percent.  The  total  error  in  the  relative  cryogenic  measure- 
ments is  estimated  to  be  no  greater  than  + 4 percent. 

IV.  THE  STRUCTURE  OF  FUSED  SILICA 

Extensive  x-ray  diffraction  studies  of  fused  silica  conducted 
primarily  by  h'arren  (1937)  have  confirmed  the  structural  arrangement 
suggested  by  Zacharicsen  (1932).  These  studies  reveal  that  the  basic 
structure  of  fused  silica  is  that  of  a silicon  atom  tetrahedral lv 
surrounded  by  four  oxygen  atoms.  The  tetrahedra  are  then  linked  together 
through  the  Si -0-Si  bonds  in  a random  network  arrangement.  Breaks  in 
the  structure  occur  every  4th  to  20th  tetrahedron  (Primak,  1975)  and  are 

o 

often  an  impurity  site.  The  average  Si-0  distance  is  about  1.6  A and  the 

o 

average  Si-Si  distance  is  about  3.0  A. 

The  Si-0  bonding  is  neither  covalent  nor  electrovalent  but  a 
combination  of  both  (Simon,  1958).  The  electrovalent  character  of  the 
Si-0  bonds  are  responsible  for  rapid  reforming  of  broken  bonds  but  the 
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covalent  character  result,  in  a structure  which  is  not  rcadilv 
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rearranged.  The  covalent  character  also  leads  to  a normal  Si hojul 
angle  of  approximately  140°. 


V.  DISCUSSION  OF  Till;  ROOM  TEMPERATURE 
EXPERIMENTAL  RESULTS 


An  examination  of  Table  V. 2,  page  72,  reveals  that  the  second-order 
elastic  constants  of  the  fused  silica  samples  with  high  Oil  content 
(Suprasil  1 and  Suprasil  2)  are  lower  than  the  samples  with  low  Oil  con- 
tent (Suprasil  Ivl  and  Suprasil  W2).  These  results  are  consistent  with 
those  of  Iletherington  and  Jack  (1962)  who  noted  also  that  an  increase 
in  OH  content  caused  an  increase  in  mass  density  fluctuations 
(inhomogeneity)  as  much  as  3 parts  in  10A  in  their  experiments.  Primak 
(1975)  suggested  that  in  similar  measurements  by  Eraser  (1968)  in  tie 
ultrasonic  range  of  frequencies  the  results  reflect  the  effect  of  no- - 
wo'-k  cleavages  induced  by  the  OH  impurities. 

Table  V.2  shows  that  8 is  negative  for  all  the  fused  silica  samples. 
This  means  that  An  is  also  negative  and  (from  Eq.  (V.4)]  that  C...  is 

t.  ill 

positive.  Bains  and  Breazeale  (1975)  have  shown,  using  a phase 
sensitive  detector,  that  the  harmonics  of  fused  silica  are  generated 
out  of  phase  with  those  generated  by  the  same  ultrasonic  wave  in  copper. 
They  conclude  that  the  positive  of  fused  silica  means  that  the 

solid  becomes  softer  upon  compression  in  contrast  to  copper  which 
becomes  stiffer  upon  compression. 

According  to  Table  V.2  Suprasil  1 has  a substantially  higher  value 
of  Cju  than  any  of  the  other  samples  used  in  these  experiments.  It  is 
tempting  to  infer  that  this  is  because  of  the  high  OH  content  of  the 
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sample  (1200  ppm  compared  to  5 ppm  for  Suprasil  IV 1 and  Suprasil  W21, 
but  Suprasil  2 also  has  high  Oil  content  (1200  ppm)  and  gives  a value 
of  Cju  which  lies  between  those  of  Suprasil  W1  and  Suprasil  W2.  A 
check  for  residual  stresses  in  all  the  samples  with  crossed  polaroid 
sheets  revealed  that  all  of  the  samples  were  stress  free.  Hence,  the 
situation  cannot  be  explained  by  residual  stresses  in  a particular 
sample.  A possible  explanation  may  be  that  the  Suprasil  1 sample  is 
strongly  homogeneous  in  all  directions  whereas  the  Suprasil  2 sample 
is  strongly  homogeneous  only  in  the  direction  perpendicular  to  the 
sample  faces  (see  Chapter  IV,  page  56).  This  means  that  mass  density 
variations  may  exist  across  the  face  of  the  ultrasonic  beam  large 
enough  to  scatter  large  numbers  of  phonons.  Since  the  amount  of  phonon 
scattering  would  be  frequency  dependent,  fundamental  and  second  harmonic 
phonons  would  be  scattered  in  different  amounts  and  could  be  such  as  to 
cause  a decrease  in  the  value  of  for  Suprasil  2 to  that  in 

Table  V.2,  page  72.  If  this  is  the  case  then  from  Table  V.2  the  values 
of  Cjjj  for  the  low  Oil  content  samples  are  lower  than  the  values 

for  the  high  OH  content  samples.  This  is  in  contrast  to  the  values 
for  the  low  OH  content  samples  being  higher  than  the  high  OH  content 
samples.  Thus,  the  network  cleavages  introduced  by  the  OH  impurities 
which  may  be  the  cause  of  lower  values  may  also  be  the  cause  of  an 
increase  in  the  values. 

VI.  DISCUSSION  OF  TEMPERATURE  DEPENDENT 
RESULTS  FOR  Cm 

Figures  V.13-V.16,  pages  81-84,  show  a significant  amount  of 
scatter  in  the  data  points  of  as  a function  of  temperature. 


ss 


Considering  that  the  estimate  of  random  error  due  to  temperature  change* 
in  the  transducer  bonds  was  ♦ 3 percent  and  that  the  standard  error  was 
no  more  than  t 1 percent,  the  considerable  amount  of  scatter  would  seem 
to  indicate  structural  changes  in  the  samples  as  a function  of 
temperature.  This  is  consistent  with  the  suggestions  of  Primak  (197j) 
who  asserts  that  variations  do  occur  in  the  structure  of  fused  silica 
as  the  temperature  is  changed.  This  is  in  largo  part  due  to  the 
asymmetrical  nature  of  the  void  surrounding  the  oxygen  atom  which  is  a 
reflection  of  the  random  nature  of  the  tetrahedral  network.  As  the 
temperature  is  lowered  evidence  exists  (Primak,  1975)  that  the  Si-O-Si 
bond  angle  changes  due  to  the  asymmetrical  vibrations  of  the  oxygen 
atoms.  Thus,  the  tetrahedral  network  of  fused  silica  below  the  quenching 
temperature  (the  temperature  at  which  the  melt  turns  into  a solid)  is 
in  a state  of  stress.  The  resulting  states  of  strain  arc  more  accentuated 
for  the  weak  impurity  bonds  and  this  makes  them  more  efficient  electron 
traps  than  the  strong  Si-0  bonds.  This  means  that  the  effect  of 
temperature  changes  is  to  cause  variations  especially  in  .inharmonic 
dependent  effects  since  it  is  these  effects  that  are  most  sensitive  to 
variations  in  structural  symmetry  and  strain  perturbations  introduced  by 
impurity  sites. 

Although  the  systematic  and  random  errors  in  the  cryogmic 
measurements  are  too  large  to  attach  a great  amount  of  significance  to 
individual  data  points,  trends  in  the  versus  temperature  data  arc 

apparent.  The  C.^  values  of  Suprasil  W1 , Suprasil  W2,  and  Suprasil  2 
all  show  a general  decrease  of  with  decreasing  temperature.  The 

Cjjj  value  of  Suprasil  1 on  the  other  hand  certainly  does  not  generally 
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decrease  with  decreasing  temperature  and  ma\  even  exhibit  .1  slight 
increase  with  decreasing  temperature.  V-»  with  the  room  temperature 
measurements,  Suprasil  1 (with  high  >1!  content  seems  to  be  an  anomaly. 
Again,  however,  Suprasil  2 (also  with  high  OH  content)  seems  to  show 
the  general  (maybe  even  more  pronounced)  versus  temperature 

behavior  as  Suprasil  h]  and  Suprasil  2 (both  with  low  OH  content), 
hollowing  the  explanation  given  for  the  room  temperature  measurements, 
it  is  conceivable  that  the  lack  of  strong  homogeneity  m the  Suprasil  2 
sample  in  all  spatial  directions  gives  rise  to  large  density  fluctuations 
across  the  face  of  the  ultrasonic  beam.  The  change  of  temperature  may 
possibly  exaggerate  these  density  fluctuations  thereby  giving  rise  to  the 
observed  temperature  dependence  of  for  Suprasil  2.  The  pronounced 

decrease  of  the  C ^ , value  of  Suprasil  2 (Fig.V.lo,  page  84)  with 
decreasing  temperature  in  the  range  147  °K  to  88  °K  indicates  that  a 
strong  structural  change  (and  consequent  density  changes)  may  have  taken 
place. 

VI 1.  DISCUSSION  01  Till  Tl  MPTRATURi:  Dl: PIMM' NT  R1 'SUITS  I OR 
B AND  ITS  IMPLICATIONS  1 OR  Mil.  l ALCUl  \ HON 
OF  CrOnF.  ISLN'S  CAi  IMA 

The  plots  of  8 versus  temperature  (Figures  V.9-V.12,  pages  77-80) 
indicate  that  is  a much  weaker  function  of  temperature  for  all 

A 

the  samples  used  in  these  experiments  than  is  C^j.  Since  the  -y^tP^ ,N) 
component  of  the  strain  generalized  GrLineiscn  parameter  for  fused  silica 
is  directly  proportional  to  8 it  is  also  weakly  dependent  upon  temperature 
for  these  samples.  This  implies  that  any  significant  temperature  variation 
in  the  vo  1 time  generalized  Grune i sen  parameters  for  fused  silica 
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| IU| « i 1 1 1 . 1 a j | must  come  from  significant  variations  in  the  clastic 
constants  other  than  and 

It  was  mentioned  previously  that  the  Brunei  sen  parameter  > depends 
on  temperature  directly  through  the  modal  specific  heats  and  indirectly 
through  the  combination  of  elastic  constants  in  terms  of  which  the 
volume  generalized  Gruneisen  parameters  may  he  expressed  (Eqs.  (Ill *22, 

1 1 1. 43 1 1 . Figure  V.  17  shows  the  experimentally  determined 
Gruneisen  parameter  y as  a function  of  temperature  as  measured  In  Khitc 
(1964)  for  fused  silica  (Spcct ros i 1 ) . It  is  seen  that  y becomes  negative 
between  .2  and  .3  of  the  Debye  temperature  (495  °K  for  Spectrosil)  and 
continues  to  become  more  negative  with  decreasing  temperature.  As  seen 
from  Figure  v. 17  this  is  in  contrast  to  germanium  for  which  y levels  off 
as  the  temperature  approaches  zero  and  to  quartz  which  maintains  a 
positive  Gruneisen  parameter  as  a function  of  temperature. 

It  is  apparent  from  Eq.  (T 1 1.221  that  in  order  to  obtain  a negative 
value  of  the  Gruneisen  parameter  at  some  temperature  at  least  one  of  the 
vo 1 umc  generalized  Gruneisen  parameters  must  he  negative  at  that  temperature 

A 

since  the  modal  heat  capacities  C(p,N)  are  always  positive.  Barron  (1957) 
and  Blackman  (1958)  have  suggested  theoretical  models  which  for  some  vo 1 umc 
generalized  Gruneisen  parameters  associated  with  certain  transverse 
vibrational  modes  of  that  model  do  become  negative.  These  transverse 
modes  necessarily  involve  elastic  constants  other  than  C and  C^. 

IVhitc  and  Birch  (1965)  suggested  that  the  transverse  vibrations 

of  fused  silica  are  associated  with  the  oxygen  atoms  of  the  tetrahedral 

network. 


According  to  l.q.  (Ill 


I the  calculation  of  the  Brunei  sen 


parameter  depends  strongly  on  the  modal  heat  capacities  C(p,N)  which 
serve  as  '•weighting”  functions  of  the  y(p,N).  further,  the  denominator 
of  the  expression  is  the  Debye  heat  capacity,  Experimental Iv  determined 
measurements  of  the  heat  capacities  for  fused  silica  at  low  temperature 
reveal  that  the  experimental  values  are  larger  than  the  calculated 
Debye  heat  capacity  by  a factor  of  r>  at  certain  temperatures. 

Ilubacher,  Leadbetter,  Morrison,  and  Stoicheff  (1959)  suggested  from 
Raman  and  Brillouin  spectra  experiments  that  the  "excess'*  heat  capacity 
of  fused  silica  is  contributed  by  optical  modes  of  very  low  frequencies 
and  unknown  origin.  They  showed  that  t lie  addition  of  three  1 instein 
heat  capacity  terms  of  low  optical  mode  frequencies  to  the  calculated 
Debye  heat  capacity  gives  good  agreement  with  experiment. 

Clark  and  Strakna  (1962)  suggested  two  equilibrium  positions  for 
the  oxygen  atoms  (double  potential  well)  to  explain  the  "excess"  heat 
capacity  but  their  suggestion  fit  the  experimental  data  above  2 °K  less 
wel  than  the  optical  mode  model  of  Flubacher  ct  al.  (Leadbettcr  and 
Morrison,  1963).  However,  Arnold,  llunklinger,  Stein,  and  Dransfeld 
(1974)  have  recently  applied  the  double  potential  well  model  to  ultra- 
sonic attenuation  at  gigahertz  frequencies  in  fused  silica  below  2 °K 
and  have  found  good  qualitative  agreement  with  experiment. 

in  conclusion,  the  present  experimental  determination  of  h and 
Cjjj  as  a function  of  temperature  allows  an  important  inference  to  be 
made  regarding  the  temperature  dependence  of  the  Crime i sen  parameter 
for  fused  silica.  If  the  quasi harmonic  Debye  model  of  a solid  is 
assumed  to  hold  for  fused  silica,  then  the  relatively  weak  temperature 
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dependence  of  6 and  < jjj  indicates  that  elastic  constants  other  than 
Cj  j and  Cj  ^ must  dominate  the  strong  temperature  dependence  of  the 
experimentally  determined  Griincisc*i  parameter  for  fused  silica  at  low 
temperatures.  This  is  consistent  with  the  arguments  presented  In- 
Bar  ron  (1957),  Blackman  (1958),  White  and  Birch  (1965),  and  Primak 
(1975)  that  the  transverse  vibrations  of  the  oxygen  atoms  in  the  tetra- 
hedral network  of  fused  silica  give  rise  to  negative  values  of  the 
(irtlne i sen  parameter  at  low  temperatures.  The  elastic  constants 
associated  with  these  transverse  modes  would  he  of  the  type  Cp,  C^f 
^114'  ^456  ^or  so^^s  *n  general.  For  fused  silica  there  are  only  two 

independent  second-order  elastic  constants  (C  and  ('  , ) and  three 

11  4 ‘1 

independent  third-order  elastic  constants  (C.,.,  C and  C c ). 

Ill  1 44  4 b() 

Therefore,  onlv  C,„,  C...f  and  f would  dominate  in  the  transverse 

44  144  4b 6 

modes  for  fused  silica. 

VIII.  SUGGESTIONS  FOR  FURTHER  WORK 

The  results  of  this  study  have  raised  a number  of  questions  to  be 
investigated. 

1.  The  temperature  dependence  of  the  Griineisen  parameter  6 in 
fused  silica  cannot  he  fully  understood  until  the  elastic  constants 
other  than  anj  are  measured  as  a function  of  temperature.  The 

ultrasonic  beam  mixing  technique  of  Dunham  and  Huntington  (1970),  for 
example,  would  allow  the  complete  isolation  of  all  the  elastic  constants 
for  isotropic  solids  when  the  results  are  combined  with  the  results  of 
the  present  technique.  At  present,  the  beam  mixing  technique  has  not 
been  adapted  to  cryogenic  measurements. 
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2.  Fused  silica  is  one  of  several  solids  for  which  the  therm.!  1 
expansion  coefficients  change  sign  at  certain  temperatures.  In  the 
case  of  fused  silica  the  change  in  sign  is  associated  with  transverse 
vibrational  modes.  Water  also  exhibits  a change  in  sign  of  its  thermal 
expansion  coefficient  at  4 °C.  Since  water  is  not  known  to  support 
transverse  acoustic  waves,  it  would  be  interesting  to  investigate  the 
ultrasonic  nonlinearity  characteristics  of  water  using  the  present 
technique  around  that  temperature.  This  would  necessitate  the  develop- 
ment of  a capacitive  detector  for  liquids. 

3.  Fused  silica  is  an  anamolous  material  thermodynamically.  There 
exist  other  isotropic  materials  which  do  not  exhibit  these  anamolous 
properties  and  an  investigation  of  their  ultrasonic  nonlinearity 
properties  as  a function  of  temperature  may  yield  information  which 
would  help  to  characterize  the  difference  between  the  materials. 

4.  Adapting  the  capacitive  driver-capacitive  receiver  system  to 
the  cryogenic  apparatus  would  eliminate  any  error  introduced  by  the 
change  of  temperature  in  the  transducer  bonding  material.  It  would  also 
permit  accurate  measurements  of  ultrasonic  velocity  and  attenuation  in 
the  samples  as  functions  of  temperature. 
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APPENDIX 


THL  EFFECT  OF  DIFFRACTION  ON  THE  GATED  DOUBLL-PULS1 
SUPERPOSITION  METHOD  OF  ULTRASONIC 
VF LOC I TY  MEASUREMENTS 

The  contributions  of  diffraction  effects  to  accurate  measurements 
of  velocity  and  attenuation  of  ultrasonic  waves  is  well  known  (Seki, 
(iranato,  and  Truell,  12*^6;  Papadakis,  1972).  Correction  formulas 
(Papadakis,  1967,  1972)  for  these  diffraction  effects  have  been  derived 
for  the  ultrasonic  pulse  t ime-of-f  1 i ght  measurement  techniques  (Papadakis, 
1967;  McSkimin,  1961;  May,  1958).  The  theory  is  now  extended  to  provide 
corrections  in  velocity  measurements  utilizing  the  gated  double  pulse 
superposition  method  of  Williams  and  Lamb  (1958)  (ultrasonic  interferometer 
technique) . 

The  problem  is  that  of  a circular  piston  transducer  of  radius  a 
radiating  an  acoustic  wave  toward  a coaxial  circular  receiver  of  the  same 
dimension  and  at  a di stance  z away . If  it  is  assumed  that  the  electrical 
response  of  the  receiver  is  proportional  to  the  average  pressure  over  the 
receiver  at  a given  time,  the  problem  may  be  divided  into  two  parts: 

1.  Find  the  acoustic  pressure  field  at  an  arbitrary  point  x away 

from  the  transducer.  ' 

2.  Average  the  generated  acoustic  pressure  field  over  the  surface 
of  the  receiving  transducer. 

If  we  let  <f>  be  the  acoustic  velocity  potential,  the  acoustic 
pressure  field  may  be  found  from  the  solution  of  the  wave  equation, 
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which  satisfies  the  boundary  condition. 


^ = u()  ( constant ) , 


( A.  2) 


on  the  surface  of  the  transmitting  transducer,  c the  sound  velocity 

in  Lq.  (A.lJ,u  is  the  particle  velocity,  and  — indicates  normal  dcriva- 

u tin 

tion  in  lq.  (A. 2). 

The  acoustic  pressure  is  obtained  from 


ao 

p = po  w 


(A.  31 


where  is  the  mean  density  of  the  sound  propagation  medium.  The 
solution  to  Lq . (A.l)  may  be  written  as 


(ujt-k-  r ) 


<M  x)  = 


(A.  4) 


where  r is  the  distance  from  the  surface  element  do  ^ of  the  transmitting 
transducer  to  the  field  point  x,  k is  the  wave  number,  is  the  angular 
frequency  of  the  wave,  and  the  integral  is  performed  over  the  surface  of 
the  transmitting  transducer.  Substituting  Lq.  (A. 4)  into  fq.  (\.3)  gives 


p(x>  = -jr~  J 

o 


()  0 r sin(ut-k»r) 


. “P0U0  r cos(ujt-k*f ) , , , - 

- > ~2T~  1 r d°l  ' (,V-a 


The  acoustic  field  pressure  is  the  real  part  of  p(x). 
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Following  Scki,  (iranato,  and  True  11  (193<>)  one  may  evaluate  tlic 
real  part  of  !:.q*  (A.5)  ( in  cylindrical  polar  coordinates)  as 

p(u,z,t)  = X(u,z)  cos(u)t-kz)  ♦ Y(u,z)  sin(wt-kz)  (A.<0 

where  2 is  the  direction  perpendicular  to  the  surface  of  the  transmitting 
transducer,  u is  measured  perpendicular  to  the  z direction  at  some 
distance  z from  the  transducer,  and 


7 7 


v , i k(x~+a  ) . k(x  «-a“)  . ^ 

X(y,z)  = VQ  cos  ■ - + V sin  — - 1 (A.  ) 


..  k(x‘'+a“)  k(x“«-a“) 

H»,z)  = V0  sin  rz Vj  cos  — 


(A.S) 


am 


*„-  i i-Dno-,n 

n=0 


(A.  9) 


V,  . I (-nn{j)"2"'1 

n=0 


(A. 10) 


This  completes  part  (1)  of  the  problem. 

The  average  of  the  generated  acoustic  pressure  field  over  the 
surface  of  the  receiving  transducer  may  be  written  as 

<p(z)>  = ““  / pdon  (A. Ill 

2 °2 

where  is  the  surface  area  of  the  receiving  transducer.  Substituting 


Ft],  (A. 6)  into  (A.  11)  gives 


w 


n 


hmmbrp 


p(z)>  = — [cos (wt-kz)  / X do-  - sin((jt-kz)  / V do. 
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We  may  rewrite  (A. 12)  in  the  form 


<p(z)>  = <p>max  cos [o>t -k z-0 ( z ) ] 


(A. 13) 


where 


and 


<P 


max 


4l(/ 


«!>(  = ) = 


2 r 2 \/~> 

X do ^ ) ♦ (/  Y do2)“] 


fj 

n2 

V do. 

/ 

X do  , 

°2 

(A. 14) 


(A. 15) 


Benson  and  Kiyohara  (1974)  have  evaluated  numerically  Hqs . (A. 13) 
and  (A.  15)  as  a function  of  the  dimensionless  parameter, 


S 


z\ 

2 9 
a*- 


(A. 10) 


where  A is  the  ultrasonic  wavelength. 

lor  the  ultrasonic  interferometer  technique  we  must  consider  two 
pulse-echo  trains,  the  second  pulsc-echo  train  being  delayed  with 
respect  to  the  first  such  that  superposition  of  the  two  trains  is 
achieved.  If  attenuation  of  the  echo  signals  is  neglected,  the  total 
average  pressure  at  the  receiving  transducer  due  to  a given  pair  of 
superimposed  echoes  may  be  written  as 
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i)  ■ = ‘P  • { cos  |t.»t -k  (m*  1 ) 16-4*  ( S )|  ♦ cos  |«»t -k  (Jn*  1)  l - j*  (S  111 

1 1 1 max  m 1 1 u 1 4 ... 

l\.l  1 

where  v.  is  the  distance  between  the  transmitting  and  receiving  transducers 
(sample  length),  m and  n are  integers  representing  the  number  of  rounJ 
trips  through  the  sample  of  the  ini t ial  and  delayed  pulses  (i.e.,  the  echo 
numbers) , and 


S 

m 


~m  _ (2m+l)JU 
2 2 

a*'  a 


(A. IS) 


Thus,  we  may  write 


• p>  = 2*p>  cos  idjt-kf, (m+n+1 ) - 
1 r 1 max 


X cos 


kf.  (m-n ) 


4,(S  |-*(S  j) 


m n 


(A.  1*1) 


Nulls  in  <p>  occur  at 


k£ (m-n) 


*(SJ-<KSJ 


m n 'u+ 


= (^v-)  * 


(A. JO) 


where  q = integer.  If  we  consider  two  nulls  of  Iq.  (A.  19)  at  q ^ and  q, 
and  write 


k = 


2tt  f 


( A . 2 1 ) 


where  f is  the  ultrasonic  frequency  and  v is  the  ultrasonic  velocity, 
then  we  may  write  from  I.q.  (A. JO) 


100 


